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Abstract: We describe a program to prove the Kepler conjecture on sphere pack- 
ings. We then carry out the first step of this program. Each packing determines 
a decomposition of space into Delaunay simplices, which are grouped together into 
finite configurations called Delaunay stars. A score, which is related to the density 
of packings, is assigned to each Delaunay star. We conjecture that the score of 
every Delaunay star is at most the score of the stars in the face-centered cubic and 
hexagonal close packings. This conjecture implies the Kepler conjecture. To com- 
plete the first step of the program, we show that every Delaunay star that satisfies 
a certain regularity condition satisfies the conjecture. 



Contents: 1. Introduction, 2. The Program, 3. Quasi-regular Tetrahedra, 4. 
Quadrilaterals, 5. Restrictions, 6. Combinatorics, 7. The Method of Subdivi- 
sion, 8. Explicit Formulas for Compression, Volume, and Angle, 9. Floating-Point 
Calculations. 

Appendix. D. J. Muder's Proof of Theorem 6.1. 

Section 1. Introduction 



The Kepler conjecture asserts that no packing of spheres in three dimensions has 
density exceeding that of the face-centered cubic lattice packing. This density is 
tt/vTs ~ 0.74048. In an earlier paper [H2], we showed how to reduce the Kep- 
ler conjecture to a finite calculation. That paper also gave numerical evidence in 
support of the method and conjecture. This finite calculation is a series of optimiza- 
tion problems involving up to 53 spheres in an explicit compact region of Euclidean 
space. Computers have little difficulty in treating problems of this size numerically, 
but a naive attempt to make a thorough study of the possible arrangements of these 
spheres would quickly exhaust the world's computer resources. 

The first purpose of this paper is to describe a program designed to give a rigorous 
proof of the Kepler conjecture. A sketch of a related program appears in [H2]. 
Although the approach of [H2] is based on substantial numerical evidence, some of 
the constructions of that paper are needlessly complicated. This paper streamlines 
some of those constructions and replaces others with constructions that are more 
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amenable to rigorous methods. For this program to succeed, the original optimiza- 
tion problem must be partitioned into a series of much smaller problems that may 
be treated by current computer technology or hand calculation. 

The second purpose of this paper is to carry out the first step of the proposed 
program. A statement of the result is contained in Theorem 1 below. 

Background to another approach to this problem is found in [H3]. To add more 
detail to the proposed program, we recall some constructions from earlier papers 
[HI] , [H2] . Begin with a packing of nonoverlapping spheres of radius 1 in Euclidean 
three-space. The density of a packing is defined in [HI] . It is defined as a limit of the 
ratio of the volume of the unit balls in a large region of space to the volume of the 
large region. The density of the packing may be improved by adding spheres until 
there is no further room to do so. The resulting packing is said to be saturated. It 
has the property that no point in space has distance greater than 2 from the center 
of some sphere. 

Every saturated packing gives rise to a decomposition of space into simplices called 
the Delaunay decomposition. The vertices of each Delaunay simplex are centers 
of spheres of the packing. None of the centers of the spheres of the packing lie 
in the interior of the circumscribing sphere of any Delaunay simplex. In fact, this 
property is enough to completely determine the Delaunay decomposition except for 
certain degenerate packings. A degeneracy occurs, for instance, when two Delaunay 
simplices have the same circumscribing sphere. In practice, these degeneracies 
are important, because they occur in the face-centered cubic and hexagonal close 
packings. The paper [H2] shows how to resolve the degeneracies by taking a small 
perturbation of the packing. In general, the Delaunay decomposition will depend 
on this perturbation. Wc refer to the centers of the packing as vertices., since the 
structure of the simplicial decomposition of space will be our primary concern. For 
a proof that the Delaunay decomposition is a dissection of space into simplices, we 
refer the reader to [R] . 

The Delaunay decomposition is dual to the well-known Voronoi decomposition. If 
the vertices of the Delaunay simplices are in nondegenerate position, two vertices 
arc joined by an edge exactly when the two corresponding Voronoi cells share a 
face, three vertices form a face exactly when the three Voronoi cells share an edge, 
and four vertices form a simplex exactly when the four corresponding Voronoi cells 
share a vertex. In other words, two vertices are joined by an edge if they lie on a 
sphere that does not contain any other of the vertices, and so forth (again assuming 
the vertices to be in nondegenerate position). The collection of all simplices that 
share a given vertex is called a Delaunay star. (This is a provisional definition: it 
will be refined below.) 

Every Delaunay simplex has edges between 2 and 4 in length and, because of the 
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simplex S in this paper comes with a fixed order on its edges, 1, . . . ,6. The order 
on the edges is to be arranged so that the first, second, and third edges meet at 
a vertex. We may also assume that the edges numbered i and i + 3 are opposite 
edges for i — 1,2, 3. We define S{yi, . . . , i/q) to be the (ordered) simplex whose ith 
edge has length j/j. If 5 is a Delaunay simplex in a fixed Delaunay star, then it 
has a distinguished vertex, the vertex common to all simplices in the star. In this 
situation, we assume that the edges are numbered so that the first, second, and 
third edges meet at the distinguished vertex. 

A function, known as the compression T{S), is defined on the space of all Delaunay 
simplices. Let Soct = (-Stt + 12 arccos(l/v/3))/v^ ^ 0.720903 be the density of a 
regular octahedron with edges of length 2. That is, place a unit ball at each vertex 
of the octahedron, and let doct be the ratio of the volume of the part of the balls 
in the octahedron to the volume of the octahedron. Let 5" be a Delaunay simplex. 
Let B be the union of four unit balls placed at each of the vertices of S. Define the 
compression as 

T{S) = -6octvol{S) + Yol{S n B). 

We extend the definition of compression to Delaunay stars D* by setting r{D*) = 
^r(5'), with the sum running over all the Delaunay simplices in the star. 

In this and subsequent work, we single out for special treatment the edges of length 
between 2 and 2.51. The constant 2.51 was determined experimentally to have a 
number of desirable properties. This constant will appear throughout the paper. 
We will call vertices that come within 2.51 of each other close neighbors. 

We say that the convex hull of four vertices is a quasi-regular tetrahedron (or sim- 
ply a tetrahedron) if all four vertices are close neighbors of one another. Suppose 
that we have a configuration of six vertices in bijection with the vertices of an 
octahedron with the property that two vertices are close neighbors if and only if 
the corresponding vertices of the octahedron are adjacent. Suppose further that 
exactly one of the three diagonals has length at most 2\/2. In this case we call the 
convex hull of the six vertices a quasi-regular octahedron (or simply an octahedron) . 

The compatibility of quasi-regular tetrahedra and octahedra with the Delaunay 
decomposition is established in Section 3. We think of Euclidean space as the 
union of quasi-regular tetrahedra, octahedra, and various less interesting Delaunay 
simplices. From now on, a Delaunay star is to be the collection of all quasi-regular 
tetrahedra, octahedra, and Delaunay simplices that share a common vertex v. This 
collection of Delaunay simplices and quasi-regular solids is often, but not always, 
the same as the objects called Delaunay stars in [H2]. We warn the reader of this 
shift in terminology. 
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It is convenient to measure the compression in multiples of the compression of 
the regular simplex of edge length 2. We define a point (abbreviated pt) to be 
r(5(2,2,2,2,2,2)). We havept= IItt/S - 12 arccos(l/^/3) ^ 0.0553736. 

One of the main purposes of this paper and its sequel is to replace the compression 
by a function (called the score) that has better properties than the compression. 
Further details on the definition of score will appear in Section 2. We are now able 
to state the main theorem of this paper. 

Theorem 1. If a Delaunay star is composed entirely of quasi- regular tetrahedra, 
then its score is less than 8 pt. 

The idea of the proof is the following. Consider the unit sphere whose center is 
the center of the Delaunay star D* . The intersection of a simplex in D* with this 
unit sphere is a spherical triangle. For example, a regular tetrahedron with edges 
of length 2 gives a triangle on the unit sphere of arc length tt/S. The star D* gives 
a triangulation of the unit sphere. The restriction on the lengths of the edges of a 
quasi-regular tetrahedron constrains the triangles in the triangulation. We classify 
all triangulations that potentially come from a star scoring more than 8 pt. Section 
5 develops a long list of properties that must be satisfied by the triangulation of a 
high-scoring star. 

It is then necessary to classify all the triangulations that possess the properties 
on this list. The original classification was carried out by a computer program, 
which generated all potential triangulations and checked them against the list. D. 
J. Muder has made a significant improvement in the argument by giving a direct, 
computer-free classification. His result appears in the appendix. 

As it turns out, there is only one triangulation that satisfies all of the properties 
on the list. Section 7 proves that Delaunay stars with this triangulation score less 
than 8pt. This will complete the main thread of the argument. 

There are a number of estimates in this paper that are established by computer. 
These estimates are used throughout the paper, even though their proofs are not 
discussed until Sections 8 and 9. These sections may be viewed as a series of 
technical appendices giving explicit formulas for the compression, dihedral angles, 
solid angles, volumes, and other quantities that must be estimated. The final section 
states the inequalities and gives details about the computerized verification. There 
is no vicious circle here: the results of Sections 8 and 9 do not rely on anything 
from Sections 2-7. 

There are several functions of a Delaunay simplex that will be used throughout this 
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defined to be the dihedral angle of the simplex S along the first edge (with respect 
to the fixed order on the edges of S). Set dih^in = dih(5(2, 2, 2.51, 2, 2, 2.51)) ^ 
0.8639; dih^^ax = dih(5(2.51, 2, 2, 2.51, 2, 2)) = arccos(-29003/96999) ^ 1.874444. 
We will see that dihmin and dih^ax are lower and upper bounds on the dihedral 
angles of quasi-regular tetrahedra. The solid angle (measured in steradians) at the 
vertex joining the first, second, and third edges is denoted so^S"). The intersection 
of S with the ball of unit radius centered at this vertex has volume sol(S')/3 (see 
[Hl,2.1]). For example, sol(5(2, 2, 2, 2, 2, 2)) fti 0.55. Let rad(^) be the circumradius 
of the simplex S. In Section 2, we will define two other functions: vor(5'), which 
is related to the volume of Voronoi cells, and the score cr{S). Finally, let r]{a,b,c) 
denote the circumradius of a triangle with edges a, b, c. Explicit formulas for all 
these functions appear in Section 8. 

Section 2. The Program 

By proving Theorem 1, the main purpose of this paper will be achieved. Neverthe- 
less, it might be helpful to give a series of comments about how Theorem 1 may be 
viewed as the solution to the first of a handful of optimization problems that would 
collectively provide a solution to the Kepler conjecture. 

We begin with some notation and terminology. We fix a Delaunay star D* about a 
vertex vq, which we take to be the origin, and we consider the unit sphere at vq. Let 
vi and V2 be vertices of D* such that fo, fi, and V2 are all close neighbors of one 
another. We take the radial projections pi of Vi to the unit sphere with center at the 
origin and connect the points pi and p2 by a geodesic arc on the sphere. We mark 
all such arcs on the unit sphere. Lemma 3.10 will show that the arcs meet only 
at their endpoints. The closures of the connected components of the complement 
of these arcs are regions on the unit sphere, called the standard regions. We may 
remove the arcs that do not bound one of the regions. The resulting system of edges 
and regions will be referred to as the standard decomposition of the unit sphere. 

Let C be the cone with vertex vq over one of the standard regions. The collection 
of the Delaunay simplices, quasi-regular tetrahedra, and quasi-regular octahedra of 
D* in C (together with the distinguished vertex vq) will be called a standard cluster. 
Each Delaimay simplex in D* belongs to a unique standard cluster. Each triangle 
in the standard decomposition of the unit sphere is associated with a unique quasi- 
regular tetrahedron, and each tetrahedron determines a triangle in the standard 
decomposition (Lemma 3.7). We may identify quasi-regular tetrahedra with clusters 
over triangular regions. 

We assign a score to each standard cluster in [H4,3]. In this section we define 
the score of a quasi-regular tetrahedron and describe the properties that the score 
should have in general. 
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Let 5 be a quasi-regular tetrahedron. It is a standard cluster in a Delaunay star 
with center vq. If the circumradius of S is at most 1.41, then we define the score 
to be T{S). 

If the circumradius is greater than 1.41, then embed the simplex S in Euclidean 
three-space. Partition Euclidean space into four infinite regions (infinite Voronoi 
cells) by associating with each vertex of S the points of space closest to that vertex. 
By intersecting S with each of the four regions, we partition S into four pieces ^o. 
Si, S2, and S3, corresponding to its four vertices vq, vi, V2, and v^. Let solj be the 
solid angle at the vertex Vi of the simplex. The expression — 45octvol(/S'o) +4 solo /3 
is an analytic function of the lengths of the edges for simplices S that contain 
their circumcenters. (Explicit formulas appear in Section 8.) This function may 
be analytically continued to a function of the lengths of the edges for simplices 5" 
that do not necessarily contain their circumcenters. Let vor(5') be defined as the 
analytic continuation of — 45octvol(5'o) + 4solo /3. 

In this case, define the score of S to be vor(S'). Write cr{S) for the score of a 
quasi-regular tetrahedron. In summary, the score is 

f r(5'), if the circumradius of S is at most 1.41, 
a(S) = \ ^ \ 

{ vor(S'), otherwise. 



The quasi-regular tetrahedron S appears in four Delaunay stars. In the other three 
Delaunay stars, the distinguished vertices will be vi, V2, and ^3, so that 5", viewed as 
a standard cluster in the other Delaunay stars, will have scores —Soct'vol{Si) +soli /3 
(or their analytic continuations), for i — 1,2, 3. By definition, 

3 

^(-5) = J2{-Soctyol{Si) + soli /3). 

i=0 

The sum of the scores of S, for each of the four vertices of S, is 4r(S'). This is the 
same total that is obtained by summing the compression of S at each of its vertices. 
This is the property we need to relate the score to the density of the packing. It 
means that although the score reapportions the compression among neighboring 
Delaunay stars, the average of the compression over a large region of space equals 
the average of the score over the same region, up to a negligible boundary term. 

The analytic continuation in the definition of vor(jS') has the following geometric 
interpretation. If the circumcenter of a quasi- regular tetrahedron S is not contained 
in 5, a small tip of the infinite Voronoi cell at vq (or some other vertex) will 
protrude through the opposite face of the Delaunay simplex. The volume of this 
small protruding tip is not counted in — 45octVol(S'o) -|-4solo /3, but it is counted in 
the analytic continuation. The analytic continuations of the scores of S for each of 
the other three vertices acquire a term representing the negative volume of a part 
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of the tip. The three parts together constitute the entire tip, so that the negative 
volumes exactly offset the volume of the tip, and the sum of the four scores of S is 
4r(S'). Details appear in [H4]. 

The general definition of the score will have similar properties. To each standard 
cluster of a Delaunay star D* a score will be assigned. The rough idea is to let the 
score of a simplex in a cluster be the compression r{S) if the circumradius of every 
face of S is small, and otherwise to let the score be defined by Voronoi cells (in a 
way that generalizes the definition for quasi- regular tetrahedra). 

The score cr(Z)*) of a Delaunay star is defined as the sum of the scores of its standard 
clusters. The score has the following properties [H4,3.1 and 3.5]. 

1. The score of a standard cluster depends only on the cluster, and not on the way 
it sits in a Delaunay star or in the Delaunay decomposition of space. 

2. The Delaunay stars of the face-centered cubic and hexagonal close packings score 
exactly 8 pt. 

3. The score is asymptotic to the compression over large regions of space. We 
make this more precise. Let A denote the vertices of a saturated packing. Let A^v 
denote the vertices inside the ball of radius N. (Fix a center.) Let D*{v) denote 
the Delaunay star at v E A. Then the score satisfies (in Landau's notation) 

Ajv Aiv 

Lemma 2.1. If the score of every Delaunay star in a saturated packing is at most 
s < 167r/3, then the density of the packing is at most 167r5oct/(167r — 3s). If the 
score of every Delaunay star in a packing is at most 8pt, then the density of the 
packing is at most tt / \/TS. 

Proof: The second claim is the special case s = 8pt. The proof relies on property 
3. The number of vertices such that D*{v) meets the boundary of the ball Bn of 
radius N has order 0{N'^). Since the Delaunay stars give a four- fold cover of IR^, 
we have 

4-7r 

4{-Soctyol{BN) + |Aiv| y) = J2i-Soctyol{D*{v)) +4sol{D%v))/3) + 0{N^) 

Ajv 

= J2r{D%v)) + 0{N') 

Ajv 

= Y,a{D*{v)) + 0{N') 

Ajv 

< s\An\+0{N^). 
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Rearranging, we get 



3vol(Sjv) ~ (1 - 3s/167r) vol(Sjv) ' 

In the limit, the left-hand side is the density and the right-hand side is the bound. 
Similar arguments can be found in [HI] and [H2]. I 



The following conjecture is fundamental. By the lemma, this conjecture implies 
the Kepler conjecture. The lemma also shows that weaker bounds than 8 pt on the 
score might be used to give new upper bounds on the density of sphere packings. 



Conjecture 2.2. The score of every Delaunay star is at most 8pt. 

The basic philosophy behind the approach of this paper is that quasi-regular tetra- 
hedra are the only clusters that give a positive score, standard clusters over quadri- 
lateral regions should be the only other clusters that may give a score of zero, and 
every other standard cluster should give a negative score. Moreover, we will prove 
that no quasi-regular tetrahedron gives more than 1 pt. 

Thus, heuristically, we try to obtain a high score by including as many triangular 
regions as possible. If we allow any other shape, preference should be given to 
quadrilaterals. Any other shape of region should be avoided if possible. If these 
other regions occur, they should be accompanied by additional triangular regions 
to compensate for the negative score of the region. We will see later in this paper 
that even triangular regions tend to give a low score unless they are arranged to 
give five triangles around each vertex. 

The main steps in a proof of the Kepler conjecture are 

1. A proof that even if all regions are triangular the total score is less than 8pt 

2. A proof that standard clusters in regions of more than three sides score at most 

Opt 

3. A proof that if all of the standard regions are triangles or quadrilaterals, then 

the total score is less than 8pt (excluding the case of pentagonal prisms) 

4- A proof that if some standard region has more than four sides, then the star 
scores less than 8pt 

5. A proof that pentagonal prisms score less than 8 pt 
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The division of the problem into these steps is quite arbitrary. They were originally 
intended to be steps of roughly equal magnitude, although is has turned out that 
a construction in [H4] has made the second step substantially easier than the long 
calculations of the third step. 

This paper carries out the first step. The second step of the program is also complete 
[H4]. Partial results are known for the third step [H5]. In the fourth step, it will 
be necessary to argue that these regions take up too much space, give too little in 
return, and have such strongly incompatible shapes that they cannot be part of a 
winning strategy. 

To make step 5 precise, we define pentagonal prisms to be Delaunay stars whose 
standard decomposition has ten triangles and five quadrilaterals, with the five 
quadrilaterals in a band around the equator, capped on both ends by five triangles 
(Diagram 2.3.) The conjecture in this section asserts, in particular, that pentagonal 
prisms, which created such difficulties in [H2], score less than 8pt. The final step 
has been separated from the third step, because the estimates are expected to be 
more delicate for pentagonal prisms than for a general Delaunay star in the third 
step. 



Diagram 2.3 




One of the main shortcomings of the compression is that pentagonal prisms have 
compression greater than 8pt (see [H2]). Numerical evidence suggests that the 
upper bound on the compression is attained by a pentagonal prism, denoted D*^^^ 
in [H2], at about 8.156 pt, and this means that the link between the compression and 
the Kepler conjecture is indirect. The score appears to correct this shortcoming. 

What evidence is there for the conjecture and program? They come as the result of 
extensive computer experimentation. I have checked the conjecture against much 
of the data obtained in the numerical studies of [H2,9.3]. The data suggest that the 
score tends to give a dramatic improvement over the compression, often improving 
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the bound by more than a point. The score of the particularly troublesome pen- 
tagonal prism -Dpp drops safely under 8 pt. I have checked a broad assortment of 
other pentagonal prisms and have found them all to score less than 8 pt. 

The second step shows that no serious pathologies can arise. The only way to form 
a Delaunay star with a positive score is by arranging a number of quasi-regular 
tetrahedra around a vertex (together with other standard clusters that can only 
lower the score). There must be at least eight tetrahedra to score 8pt, and if there 
are any distortions in these tetrahedra, there must be at least nine. However, as 
this paper shows, too many quasi-regular tetrahedra in any star are also harmful. 
Future papers will impose additional limits on the structure of the optimal Delaunay 
star. 

Section 3. Quasi-Regular Tetrahedra 

This section studies the compatibility of the Delaunay simplices and the quasi- 
regular solids. Fix three vertices vi, V2, and vs that are close neighbors to one 
another. Let T be the triangle with vertices Vi. It does not follow that T is the face 
of a Delaunay simplex. However, as we will see, when T is not the face of a simplex, 
the arrangement of the surrounding simplices is almost completely determined. 

If T is not the face of a Delaunay simplex, then we will show that there are two 
additional vertices vq and v'q, where vq and v'q are close neighbors to vi, V2, and 
V3. This means that there are two quasi- regular tetrahedra 5*1 and S2 with ver- 
tices (vo, ■i'l, '^2, 'i^s) and {vq,vi,V2,vs), respectively, that have the common face T 
(see Diagram 3.1. a). We will see that U ^2 is the union of three Delaunay 
simplices with vertices {vo,Vq,vi,V2), {vq,Vq,V2,vs), and (tJo, t'o, "i^s, "i^i) in Diagram 
S.l.b. This section establishes that this is the only situation in which quasi- regular 
tetrahedra are not Delaunay simplices: they must come in pairs and their union 
must be three Delaunay simplices joined along a common edge. The decomposition 
of this paper is obtained by taking each such triple of Delaunay simplices (S.l.b) 
and replacing the triple by a pair of quasi- regular tetrahedra (3.1. a). 
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Diagram 3.1 
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From the dual perspective of Voronoi cells, the Voronoi cell at vq (or v'q) would 
have a small tip protruding from Si through T, if the vertex v'q were not present. 
The vertex v'q slices off this protruding tip so that the Voronoi cells at vq and v'q 
have a small face in common. 

Lemma 3.2. Suppose that the circumradius of the triangle T is less than -\/2. 
Then T is the face of a Delaunay simplex. 

Proof: Let r < \/2 be the radius of the circle that circumscribes T, and let c be 
the center of the circle. The sphere of radius r at c does not contain any vertices 
of D* other than vi, ^2, and ^3. By the definition of the Delaunay decomposition 
(as described in the introduction) , this implies that T is the face of a simplex. I 



Remark 3.3. We have several constraints on the edge lengths, if T is not a face of 
a Delaunay simplex. Consider the circumradius ?7(a, 6, c) of a triangle whose edges 
have lengths o, 6, c between 2 and 2.51. Since 2.51^ < 2^+2^, we see that the triangle 
is acute, so that ?7(a, 6, c) is monotonically increasing in a, 6, and c. This gives simple 
estimates relating the circumradius to a, 6, and c. The circumradius is at most 
?7(2.51, 2.51, 2.51) = 2.51/^3 ^ 1.449. If the circumradius is at least \/2 ^ 1.41421, 
then a, 6, and c are greater than 2.3 (?7(2.3, 2.51, 2.51) ^ 1.41191 < v^). Under the 
same hypothesis, two of a, 6, and c are greater than 2.41 (77(2. 41, 2.41, 2.51) < -\/2). 
Finally, at least one edge has length greater than 2.44 (7/(2.44, 2.44, 2.44) < \f2). 

Lemma 3.4. Let T be the triangle with vertices vi. Assume that the vertices Vi 
are close neighbors of one another. Suppose there is a vertex vq that lies closer 
to the circumcenter of T than the vertices of T do. Then the vertex vq satisfies 
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2 < |fo — f i I < 2.15, for i = 1,2, 3. In particular, the convex hull of vq, . . . , is a 
quasi-regular tetrahedron S with face T. 

Another way of stating the hypothesis on the circmncenter is to say that the plane 
of T separates vq from the circumcenter of S. Because of the constraints on the edge 
lengths in Remark 3.3, the three other faces of S are faces of Delaunay simplices. 

Proof: We defer the proof to Section 8.2.5. I 

Lemma 3.5. Let v, vi, V2, V3, and V4 be distinct vertices with pairwise distances at 
least 2. Suppose that the pairs {vi,Vj) are close neighbors for ^ {1)4}. Then 

V does not lies in the convex hull of {vi, V2iV^,V4). 

Proof: For a contradiction, suppose v lies in the convex hull. Since |f — fi| > 2 and 
for 7^ {I74}, Vi and Vj are close neighbors, the angle formed by Vi and Vj at 

the vertex v is at most ^0 = arccos(l — 2.51^/8) ~ 1.357. Each such pair of vertices 
gives a geodesic arc of length at most 6*0 radians on the unit sphere centered at 
V. We obtain in this way a triangulation of the unit sphere by four triangles, two 
with edges of length at most 6q = 2 arcsin(2.51/4) radians, and two others that fit 
together to form a quadrilateral with edges of at most ^0 radians. By the spherical 
law of cosines, the area formula for a spherical triangle, and [H2,6.1], each of the 
first two triangles has area at most 3 arccos(cos 9o/{l + cos ^0)) — tt ft! 1.04. By the 
same lemma, the quadrilateral has at most the of a regular quadrilateral of 
side or about 2.8. Since the combined area of the two triangles and quadrilateral 
is less than 4%, they cannot give the desired triangulation. To see that v cannot lie 
on the boundary, it is enough to check that a triangle having two edges of lengths 
between 2 and 2.51 cannot contain a point that has distance 2 or more from each 
vertex. We leave this as an exercise. I 



Corollary 3.5.1 No vertex of the packing is ever an interior point of a quasi- regular 
tetrahedron or octahedron. 

Proof: The corollary is immediate for a tetrahedron. For an octahedron, draw the 
distinguished diagonal and apply the lemma to each of the four resulting simplices.l 

Let T be a triangle of circumradius between \/2 and 2.51 /y/S, with edges of length 
between 2 and 2.51. Consider the line through the circumcenter of T, perpendicular 
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circumcenters of the two simplices with face T formed by placing an additional 
vertex at distance two from the three vertices of of T on either side of the plane 
through T. 

Lemma 3.6. Let be a simplex formed by the vertices of T and a fourth vertex 
Vq. Suppose that the circumcenter of S lies on the segment s. Assume that v'q has 
distance at least 2 from each of the vertices of T. Then v'q has distance less than 
2.3 from each of the vertices of T. 

Proof: Let vi, V2, be the vertices of T. For a contradiction, assume that v'q 
has distance at least 2.3 from a vertex v\ of T . Lemma 3.4 shows that the plane 
through T does not separate v'q from the circumcenter of S. If the circumcenter 
lies in s (and is not separated from v'q by the plane through T), then by moving v'q 
to decrease the circumradius, the circumcenter remains in s. 

Let S be the simplex with vertices Vq, vi, vi-, and vj,- The circumcenter of S lies in 
the interior of S. We omit the proof, because it is established by methods similar 
to (but longer than) the proof of Lemma 3.4. Thus, the circumradius is increasing 
in the lengths of \v'q — f j|, for i = 1, 2, 3 (see 8.2.4). 

Let R be the circumradius of a simplex with face T and center an endpoint of 
s. We will prove that the circumradius of S is greater than i?, contrary to our 
hypothesis. Moving v'q to decrease the circumradius, we may take the distances 
to Vi to be precisely 2.3, 2, and 2. We may move vi, V2, and v^ along their fixed 
circumscribing circle until \vi — vz\ = 2.51 and \v\ — V2\ = \vi — v^\ in a way 
that does not decrease any of the distances from v'q to Vi. Repeating the previous 
step, we may retain our assumption that v'q has distances exactly 2.3, 2, and 2 
from the vertices Vi as before. We have reduced the problem to a one-dimensional 
family of tetrahedra parametrized by the radius r of the circumscribing circle of 
T. Set a;(r) = 1^1 — V3I = 1^2 — vsl- To obtain our desired contradiction, we must 
show that the circumradius R'{r) of the simplex 5'(2.3, 2, 2, 2.51, x(r), x{r)) satisfies 
R'{r) > R{r). Since both R' and R are increasing in r, for r e [V2,2.51/VS], the 
desired inequality follows if we evaluate the 200 constants 



for z = 0, . . . , 199, and check that they are all positive. (The smallest is about 



Let T be a triangle made up of three close neighbors. Suppose that T is not the face 
of a Delaunay simplex. There exists a vertex vq whose distance to the circumcenter 
of T is less than the circumradius of T. Let S be the quasi-regular tetrahedron 
formed bv vn and the vertices of T. It is not a Delaunav simolex, so there exists 





0.00005799, which occurs for i = 199.) 
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a vertex v'q that is less than the circumradius of S from the circumcenter of S. 
Let S' be the simplex formed by v'q and the vertices of T. It is not a Delaunay 
simplex either. The circumcenter of S lies in the segment s of Lemma 3.4, so the 
circumcenter of S' does too. The lengths of the edges of S and S' other than T are 
constrained by Lemmas 3.4 and 3.6. In particular, in light of Remark 3.3, the faces 
other than T oi S and S' are faces of Delaunay simplices. 

If vq and v'q lie on the same side of the plane through T, then either v'q lies in S, 
vq lies in S', or the faces of S and S' intersect nonsimplicially. None of these situ- 
ations can occur because a nondegenerate Delaunay decomposition is a Euclidean 
simplicial complex and because of Lemma 3.5. We conclude that vq and v'q lie on 
opposite sides of the plane through T. 

S U S' is bounded by Delaunay faces, so S* U S" is a union of Delaunay simplices. 
The fourth vertex of the Delaunay simplex in S U S' with face {vq, vi,V2) cannot 
be ^3 {S is not a Delaunay simplex), so it must be v'q. Similarly, {vq, v'q, V2, vs) and 
{vo,v'q,V3,vi) are Delaunay simplices. These three Delaunay simplices cannot be 
quasi- regular tetrahedra by Lemma 8.3.2. 

The assumption that T is not a face has completely determined the surrounding 
geometry: there are two quasi-regular tetrahedra S and S' along T such that SUS' 
is a union of three Delaunay simplices. 

Lemma 3.7. Let L be a union of standard regions. Suppose that the boundary of 
L consists of three edges. Then either L or its complement is a single triangle. 

For example, the interior of L cannot have the form of Diagram 3.8. Lemma 3.5 
(proof) shows that if all regions are triangles, then there are at least 12 triangles, 
so that the exterior of L cannot have the form of Diagram 3.8 either. 



Diagram 3.8 
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Proof: Replacing L by its complement if necessary, we may assume that the area of 
L is less than its complement. The triangular boundary corresponds to four vertices 
vq (the origin), vi, ^2, and v^- The close- neighbor constraints on the lengths show 
that the convex hull oivQ,... , t'3 is a quasi-regular tetrahedron. By construction 
each quasi-regular tetrahedron is a single cluster. I 



We say that a point u e is enclosed by a region on the unit sphere if the interior 
of the cone (with vertex vq) over that region contains v. For example, in Diagram 
3.9, the point v is enclosed by the given spherical triangle. 



Diagram 3.9 






V 


^^fc — — 71 / 












Vo 1^ 





The following lemma was used in Section 2 to define the standard decomposition. 

Lemma 3.10. Fix a Delaunay star D* with center vq. Draw geodesic arcs on the 
unit sphere at vq for every triple of close neighbors vq, f i, V2 (as in Section 2). The 
resulting system of arcs do not meet except at endpoints. 

Proof. Our proof will be based on the fact that a nondegenerate Delaunay decom- 
position is a Euclidean simplicial complex. Let Ti and T2 be two triangles made 
from two such triples of close neighbors. We have Ti <Z SiVJ S'^, where Si and Si 
are the Delaunay simplices with face Ti if Ti is the face of a Delaunay simplex, and 
they are the two quasi-regular tetrahedra with face Ti constructed above, otherwise. 
Since a nondegenerate Delaunay decomposition is a Euclidean simplicial complex. 
Si U S'l meets 5*2 U 5*2 simplicially. By the restrictions on the lengths of the edges in 
Lemmas 3.4 and 3.6, this forces Ti to intersect T2 simplicially. The result foUows.B 



Section 4. Quadrilaterals 

Fix a Delaunay star composed entirely of quasi-regular tetrahedra and consider the 
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by four edges of the triangulation. L will be the union of two or more triangles. 
Replacing L by its complement in the unit sphere if necessary, we assume that the 
area of L is less than that of its complement. 

Wc claim that in this context, L is the union of cither two or four triangles, as 
illustrated in Diagram 4.1. In particular, L encloses at most one vertex. If a 
diagonal to the quadrilateral L is an edge of the triangulation, the region L is 
divided into two triangles each associated with a quasi-regular tetrahedron. In 
particular, there is no enclosed vertex. (Lemma 3.7 precludes any subtriangulation 
of a triangular region.) If, however, there is a single enclosed vertex and neither 
diagonal is an edge of the triangulation, then the only possible triangulation of L 
is the one of the diagram. Proposition 4.2 completes the proof of the claim. 



Diagram 4-i 

















Proposition 4.2. A union of regions (of area less than 27r) bounded by exactly 
four edges cannot enclose two vertices of distance at most 2.51 from the origin. 

This argument is somewhat delicate: if our parameter 2.51 had been set at 2.541, 
for instance, such an arrangement would exist. First we prove a useful reduction. 

Lemma 4.3. Assume a figure exists with vectors vi, . . . , V4, v, and v' subject to 
the constraints 

2 <|t;,| < 2.51, 

2 <\vi - Vi+i \ < ki, 

2 <\v^ - Vi+2\, 

2 <\v - v'\, 
hi <\w - Vi\, 

2 <\w\ < £, for w = v,v' and z = 1, . . . , 4 (mod4) 
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where /ij, and ki are fixed constants that satisfy £ e [2.51, 2\/^], hi e [2,2\/2), 
A;j e [2, 2.51]. Let L be the quadrilateral on the unit sphere with vertices Vi/\vi\ and 
edges running between consecutive vertices. Assume that v and v' lie in the cone 
at the origin obtained by scaling L. Then another figure exists made of a (new) 
collection of vectors vi, . . . ,V4, v, and v' subject to the constraints above together 
with the additional constraints 

\vi - Vi+i\ = ki 

\vi\ = 2, for i = 1, ... ,4, 

Moreover, the quadrilateral L may be assumed to be convex. 

Proof (4.3): By rescaling v and v' , we may assume that l^;! = \v'\ = i. (Moving 
V or v' away from the origin increases its distance from the other vertices of the 
configuration.) 

The diagonals satisfy |v2— ■i'4| > 2.1 and — vs] > 2.1. Otherwise, if say V3I < 
2.1, then the faces with vertices (t'i,t'2,^3) and {vi,V4^,vs) have circumradius less 
than a/2. The edge from to v has length at most y/2, so this edge cannot intersect 
these faces by the Euclidean simplicial complex argument used before 3.7 and in 
3.10. By Lemma 3.5, v cannot lie in the convex hull of {0,vi,V3,Vi). This leaves v 
nowhere to go, and a figure with |fi — < 2.1 does not exist. 

Next, we claim that we may assume that the quadrilateral L is convex (in the sense 
that it contains the geodesic arcs between any two points in the region). To see 
this, suppose the vertex Vi lies in the cone over the convex hull of the other three 
vertices vj. Consider the plane P through the origin, Vi-i, and f^+i. The reflection 
of Vi through P is no closer to v, v' , or Vi^2 and has the same distance to the 
origin, Vi-i, and fi+i. Thus, replacing Vi with v!^ if necessary, we may assume that 
L is convex. 

Most of the remaining deformations will be described as pivots. We will fix an axis 
and rotate a vertex around a circle centered along and perpendicular to the given 
axis. If, for example, \vs — f4| < /C3, we pivot the vertex ^4 around the axis through 
and vi until |f3 — v^l — kg. It follows from the choice of axes that the distances 
from V4 to the origin and Vi are left unchanged, and it follows from the convexity of 
L that |w — f4| increases for w = V2, v^, v, and v' . Similarly, we may pivot vertices 
Vi along the axis through the origin and Vi+i until \vi — Vi-i \ — ki, for all i . 

Fix an axis through two opposite vertices (say vi and V3) and pivot another vertex 
(say V2) around the axis toward the origin. We wish to continue by picking different 
axes and pivoting until \vi\ = 2, for i = 1,2,3,4. However, this process appears to 
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vertices and the pivot toward the origin moves Vi closer to that enclosed vector. 
We must check that this situation can be avoided. 

Interchanging the roles of (^1,^3) with (t'2,'1'4) as necessary, wc continue to pivot 
until 1^1 1 = 1^3! = 2 or \v2\ = {v^l = 2 (say the former). We claim that either V2 
has distance greater than h2 from v or that pivoting V2 around the axis through vi 
and V2 moves V2 away from v. If not, we find that |fi| = |f3| —2, \v — V2\ — /i2 and 
that V lies in the cone C = C{v2) with vertex V2 spanned by the vectors from V2 to 
the origin, vi, and v^. (This relies on the convexity of the region L.) 

To complete the proof, we show that this figure made from (0, vi,V2,vs, v) cannot 
exist. Contract the edge (w, V2) as much as possible keeping the triangle (0, Vi, v^) 
fixed, subject to the constraints that v G C{v2) and — ty'l > 2, for w = v,V2 and 
w' — 0, fi, f3. This contraction gives \v — V2\ < h2 < 2\/2. 

Case 1: v lies in the plane of {0,V2,vs). This gives an impossibility: crossing 
edges {v, V2) (0, vs) of length less than 2-\/2. Similarly, v cannot lie in the plane of 

{0,Vi,V2). 

Case 2: v lies in the interior of the cone C{v2)- The contraction gives \v — w'\ = 
2 for w = v,V2 and w' = 0,vi,V3. The edge {v,V2) divides the convex hull of 
(0, fi, f2, ^3, v) into three simplices. Consider the dihedral angles of these simplices 
along this edge. The dihedral angle of the simplex {vi,V2,vs,v) is less than tt. 
The dihedral angles of the other two are less than dih(S'(2-\/2, 2, 2, 2, 2, 2)) = 7r/2. 
Hence, the dihedral angles along the diagonal cannot sum to 27r and the figure does 
not exist. 

Case 3: v lies in the plane of {v2, vi,vs). Let r be the radius of the circle in this 
plane passing through vi and f 3 obtained by intersecting the plane with a sphere 
of radius 2 at the origin. We have 2r > |f 1 — f3| > 2\/2 because otherwise we have 
the impossible situation of crossing edges (^1,^3) and {v2,v) of length less than 
2y/2. Let H be the perpendicular bisector of the segment (^1,^3). By reflecting v 
through H if necessary we may assume that v and V2 lie on the same side of H, say 
the side of ^3. Furthermore, by contracting {v,V2), we may assume without loss of 
generality that \vs — v\ = \v2 — vs\ = 2. Let /(r) = \v — as a function of r. / 
is increasing in r. The inequalities 2\/2 > h2 > \v — V2\ > f{V2) — 2\/2 give the 
desired contradiction. I 

Proof (4.2): Assume for a contradiction that v and v' are vertices enclosed by L. 
Let the center of the Delaunay star be at the origin, and let fi, . . . ,f4, indexed 
consecutively, be the four vertices of the Delaunay star that determine the extreme 
points of L. 
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We will describe a sequence of deformations of the configuration (formed by the 
vertices vi, . . . ,V4,v,v') that transform the original configuration of vertices into 
particular rigid arrangements below. We will show that these rigid arrangements 
cannot exist, and from this it will follow that the original configuration does not 
exist either. These deformations will preserve the constraints of the problem. To 
be explicit, we assume that that 2 < \w\ < 2.51, that 2 < \w — Vi\ ii w ^ Vi, that 
2 < \v — v'\, and that \vi — Vi+i\ < 2.51, for i = 1, 2, 3, 4 and w = Vi, . . . , t'4, v'. 
Here and elsewhere we take our subscripts modulo 4, so that vi — V5, and so forth. 
The deformations will also keep v and v' in the cone at the origin that is determined 
by the vertices Vj. 

We consider some deformations that increase \v — v'\. By Lemma 4.3, we may 
assume that \vi\ = 2, \v\ = \v'\ = \vi — Vi+i\ = 2.51, for i = 1,2,3,4. If, for some 
i, we have \vi — f | > 2 and \vi — v'\ > 2, then we fix Vi-i and Vi+2 and pivot Vi+i 
around the axis through the origin and f j_|_2 away from v and v' . The constraints 
\vi+i — Vi\ = \vi — Vi-i\ — 2.51 will force us to drag Vi to a new position on the 
sphere of radius 2. By making this pivot sufficiently small, we may assume that 
\vi — w\ and — w\, foTW = v, v', are greater than 2. 

The vertices v and v' cannot both have distance 2 from both and f^-i, for 
then we would have v = v' . So one of them, say v, has distance exactly 2 from 
at most one of Vi+2 and Vi-i (say ^^+2)- Thus, v may be pivoted around the axis 
through the origin and away from v' . In this way, we increase |f — f'| until 
\vi — v\ = 2 or \vi — v'\ = 2, for i = 1, 2, 3, 4. 

Suppose one of v, v' (say v) has distance 2 from Vi, Vi+i, and Vi+2- The configuration 
is completely rigid. By symmetry, the vertices Vi-i and v' must be the reflections 
of Vi+i and respectively, through the plane through 0, Vi, and In particular, 

v' has distance 2 from Vi, Vi-i, and ?Ji+2- We pick coordinates and evaluate the 
length \v — v'\. We find that \v — v'\ ~ 1.746, contrary to the hypothesis that the 
centers of the spheres of our packing are separated by distances of at least 2. Thus, 
the hypothesis that v has distance two from three other vertices is incorrect. 



(a) 



Diagram 4-4 
(b) 
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We are left with one of the configurations of Diagram 4.4. An edge is drawn in the 
diagram, when the distance between the two endpoints has the smallest possible 
value (that is, 2.51 for the four edges of the quadrilateral, and 2 for the remaining 
edges). Deform the figure of case (a) along the one remaining degree of freedom 
until |v — i;'! =2. In case (a), referring to the notation established by Diagram 4.5, 
we have a quadrilateral on the unit sphere of edges ti = 2 arcsin(2.51/4) ^ 1.357 
radians, t2 = arccos(2.51/4) radians, and ts = 2 arcsin(l/2.51) radians. The form 
of this quadrilateral is determined by the angle a, and it is clear that the angle 
P{a) is decreasing in a. We have 

9 = dihinax = dih(5(2.51, 2, 2, 2.51, 2, 2)) ^ 1.874444. 

The figure exists if and only if there exists a such that P{a) + (3{27t — 6 — a) = 
27r — 6. By symmetry, we may assume that a < tt — 9/2 k, 2.20437. The condition 
I'i^ — 'i^i I > 2 implies that 

a > arccos ((cost2 — cos t2 cos is)/ (sin ^2 sin is)) > 1.21. 

However, by monotonicity, 

(3{a) + /3(27r - 9 - a) < (3{ai) + (3{2n - 9 - ai - 0.1) <2n-9, 

for ai < a < 0.1 + a^, with = 1.21 + 0.1 z, and z = 0, 1, . . . , 9, as a direct 
calculation of the constants P{ai) + (3{2'k — 9 — — 0.1) will reveal. (The largest 
constant, which is about 2tt — 9 — 0.113, occurs for i = 0.) Hence, the figure of 
Diagram 4. 4. a does not exist. 



Diagram 4-5 
/ 27t-e-a tz \ 




t2 / 



To rule out Diagram 4.4.b, we reflect f , if necessary, to its image through the plane 
P through (0,wi,t'3), so that P separates v and v' . The vertex v can then be 
pivoted away from v' along the axis through v\ and ^3. This decreases |f |, but we 
may rescale so that \v\ = 2.51. Eventually |f — ^2] = 2 or |v — ^4] = 2. This is 
the previously considered case in which v has distance 2 from three of the vertices 
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Vi- This completes the proof that the original arrangement of two enclosed vertices 
does not exist. I 

Section 5. Restrictions 

If a Delaunay star D* is composed entirely of tetrahedra, then we obtain a trian- 
gulation of the unit sphere. As explained in Section 2, we wish to prove that no 
matter what the triangulation is, we always obtain a score less than 8pt. In this 
section we make a long list of properties that a configuration must have if it is to 
have a score of 8 pt or more. The next section and the appendix show that only one 
triangulation satisfies all of these properties. Additional arguments will show that 
this triangulation scores less than 8pt. This will complete the proof of Theorem 1. 

In this section, the term vertices refers to the vertices of the triangulation of the 
unit sphere. The edges of the triangulation give a planar graph. We adopt the 
standard terminology of graph theory to describe the triangulation. We will speak 
of the degree of a vertex, adjacent vertices, and so forth. The n triangles around 
a vertex will be referred to as an n-gon. We will also refer to the corresponding 
n tetrahedra that give the triangles of the polygon. We say that a triangulation 
contains a pattern (ai, . . . , an), for G N, if there are distinct vertices Vi of degrees 
tti that are pairwise nonadjacent, for i = 1, . . . , n. Let be the number of vertices 
in the triangulation, and let be the number of vertices of degree i. We have 
N = ZNi. 

In this section, we will start to use various inequalities from Section 9 related to 
the score. Since the score cr{S) may be either T{S) or vor(5') depending on the 
circumradius of S, there are two cases to consider for every inequality. In general, 
the inequalities for T{S) are more difiicult to establish. In the following sections, 
we will only cite the inequalities pertaining to r(5'). Section 9 shows how all the 
same inequalities hold for vor{S). 

Proposition 5.1. Consider a Delaunay star D* that is composed entirely of quasi- 
regular tetrahedra. Suppose that cr{D*) > 8pt. Then the following restrictions hold 
on the triangulation of the unit sphere given by the standard decomposition. 

1. 13<N < 15. 

2. N = + + Ne. 

3. A region bounded by three edges is either a single triangle or the complement 

of a single triangle. 

4. Two degree four vertices cannot be adjacent. 

5. A^4 < 2. 
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7. The pattern (6, 6) or (6, 4, 4) implies that N > 14. 

8. If there are two adjacent degree six vertices, and a third degree six vertex 

adjacent to neither of the first two, then N = 15 and all other vertices are 
adjacent to at least one of these three. 

9. The triangulation is made of geodesic arcs on the sphere whose radian lengths 

are between 0.8 and 1.36. 

Lemma 5.2. Consider a vertex of degree n, for some 4 < n < 7. Let 5*1, .. . , S'n 
be the tetrahedra that give the n triangles. Then ^^=ic('S'i) is less than Zn, 
where — 0.33 p^, z^ = 4.52 pt, zq = — 1.52 pt, and z-j- = — 8.9 pt. Suppose that 
n > 6. Let Si, . . . ,8^ be any four of the n tetrahedra around the vertex. Then 
Eti^iSi) < 1.5 pt. 

Proof (5.2): When n = 4, this is Lemma 9.5. When n = 5, this is Lemma 9.6. 
When n = 6 or n = 7, we have by Calculation 9.4 

n n 

^a{Si) < ^(0.378979 dih(^j) -0.410894) 

= 27r(0.378979) - 0.410894n. 

The right-hand side evaluates to about — 1.520014 and —8. 940403 pt, respectively, 
when n = 6 and n = 7. 

Assume that n > 6, and select any four Si, . . . , S4 of the n tetrahedra around the 
vertex. Let ^'5 and Sq be two other tetrahedra around the vertex. The dihedral 
angles of S5 and Sq are are at least dihmin, by Calculation 9.3. Each of the four 
triangles (associated with Si, . . . , S4) must then, on average, have an angle at most 
(27r - 2dihniin)/4 at v. By Calculation 9.4, 

4 4 

J^^C^i) < Xl(0-3'^S9^9dih(5i) - 0.410894) 

< (27r - 2dih„,in)0.378979 + 4(-0.410894) < 1.5 pt 

■ 

Proof (5.1): Let t = 2{N — 2) be the number of triangles, an even number. By 
Euler's theorem on polyhedra, 

3Ns + 2N4 + N5 + 0Ne-NT-- = 12. 

Let Si, . . . , St denote the tetrahedra of D*. Let CTj = a{Si) be the corresponding 
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47r. Often, without warning, we will rearrange the indices i so that the tetrahedra 
that give the triangles around a given vertex are numbered consecutively. When 
a vertex v of the triangulation has been fixed, we will let ctj denote the angles of 
the triangles at f , so that = 2%. The angle of a triangle is equal to the 

corresponding dihedral angle of the simplex Si. We abbreviate certain sums over n 
elements to when the context makes the indexing set clear. Throughout the 

argument, we will use Calculation 9.1, which asserts that ai < Ipt, for all i. The 
proofs will show that if a triangulation fails to have any of the properties 1-9, then 
the total score must be less than 8 pt. 

(Proof of 5.1.1): Assume that t > 28. By Calculation 9.9, 

J^(7i < ^(0.446634 soli -0.190249) < 47r(0.446634) + 1(-0.190249) 

< 47r(0.446634) + 28(-0.190249) < 8 pt. 

Suppose that t <18. By Calculation 9.8, 

J^di < J](-0.37642101soli +0.287389) = 47r(-0.37642101) + t(0.287389) 

it) (t) 

< 47r(-0.37642101) + 18(0.287389) < 8 pt. 

This proves that 12 < iV < 15. The case iV = 12 will be excluded after 5.1.2. 

(Proof of 5.1.2): By Lemma 8.3.2 and Calculation 9.3, each angle a of each triangle 
satisfies dih^in < a < dih^ax- In particular, a > 7r/4, so that each vertex has 
degree less than eight, and a < 27r/3, so that each vertex has degree greater than 
three. 

Assume for a contradiction that there is a vertex of degree seven. Consider the 
seven tetrahedra around the given vertex. By Lemma 5.2, the tetrahedra satisfy 

^C7i < -8.9 pt. 

(7) 

Suppose t < 24. For each vertex v, set = X^^cr^, the sum running over the 
tetrahedra around v. Clearly, ^(j) Ci = (^„ C)/3- Pick a vertex v that is not a 
vertex of any of the seven triangles of the heptagon. By Lemma 5.2, we see that 
Cv < 0.33 pt if V has degree four, (y < 4.52 pt if v has degree five, C„ < —1.52 pt if v 
has degree six, and (y < — 8.9 pt if v has degree seven. In particular, if v has degree 
n, then (y falls short of n points by at least 0.48 pt. Thus, 



<^^i + 5Z Pt-^-^^Pt< (-8.9+ (24- 7) - 0.48) 8;?*. 

(t) (7) (t-7) 
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Finally, we assume that t = 26 and AT = 15. If Nq = and = 1, then Euler's 
theorem gives the incompatible conditions 2N4 + N5 — I = 12 and + A^5 + 1 = 15. 
Thus, A^6 > or > 1. This gives a second A;-gon (A; = 6 or 7) around a vertex 
V. This second polygon shares at most two triangles with the original heptagon. 
This leaves at least four triangles of a second polygon exterior to the first. Thus, 
by Lemma 5.2, 

^CFi <'^ai + '^ai+ ^ pt < -8.9 pt+ 1.5 pt+ {26- 11) pt< 8 pt. 

(26) (7) (4) (t-ll) 

(Proof of 5.1.1, cont.): Assume AT = 12. We define three classes of quasi-regular 
tetrahedra. In the first class, all the edges have lengths between 2 and 2.1. In the 
second, the fourth, fifth, and sixth edges have lengths greater than 2.1. The third 
class is everything else. 

Set e = 0.001, a = -0.419351, b = 0.2856354. The following are established by 
Calculations 9.10-9.12. 

cr{S) < asol(S') + b + e, for S in the first class, 
(7{S) < asol(5') + 6, for S in the second class, 
cr{S) < asol(5') + 6 — 5e, for S in the third class. 

Consider a vertex v of degree n = 4, 5, or 6 and the surrounding tetrahedra 
Si, . . . , Sn. We claim that 

n n 

(5.1.1.1) ^(7(50 <^(asol(5i) + 6). 

i=i i=i 

This follows directly from the stated inequalities if none of these tetrahedra are in 
the first class. It is also obvious if at least one of these tetrahedra is in the third 
class, because then the inequality is violated by at most —be + (n — l)e < 0. So 
assume that all of the tetrahedra are in the first two classes with at least one in 
the first class. By the restrictions on the lengths of the edges, a tetrahedron in 
the first class cannot be adjacent to one in the second class. We conclude that the 
tetrahedra are all in the first class. 

By Lemma 5.2, Zn < n(0.904)pt. If Y.{n) soK-^i) < n(0.56176), then 

^(j{Si) <n(0.904)pt< ^(a(0.56176) + 6) < ^(a sol(5,) + 6). 

(n) (n) (n) 

So we may assume that Yl:{n) sol(5'i) > n(0. 56176). By Calculation 9.13, 

(n) (n) 

< X(asol(5i) + b) + ^(0.132365 - 0.236219(0.56176)) < X(asol(5i) + b). 
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This establishes inequality 5.1.1.1. By averaging over every vertex, we see that the 
average of the scores cr{S) is less than the average of asol{S) + b. So 

a(D*) = J2(^iS) < J^(asol(;5) + 6) =47ra + 206Ri 7.99998i?*< 8;>t 

(*) (t) 

(Proof of 5.1.3): This is Lemma 3.7. (If we had not introduced quasi- regular 
tetrahedra, then this result would no longer hold.) 

(Proof of 5.1.4): If two degree four vertices are adjacent, we then have the arrange- 
ment of Diagram 5.3. This is a quadrilateral enclosing two vertices, contrary to 
Proposition 4.2. 



Diagram 5.3 




(Proof of 5.1.5): As in the proof of 5.1.2, for each vertex v, set — J2v "^i' 
running over the tetrahedra around the vertex v. We use the estimates of that 
appear in Lemma 5.2. 

We have found that Ni = 0, if i ^ 4, 5, 6. By Euler's theorem, A^5 = 12 — 2N4 and 
Nq = N4 + N - 12. Assume that iV > 13 and that N4 > 3. Then 

^(Ti = i^Ct; < ^(0.33iV4 + 4.52(12 -2iV4)-l-52(A^4 + A^- 12)) 
< J(0.33(3) + 4.52(12 - 6) - 1.52(3 + 13 - 12)) pt < 8 pt. 

(Proof of 5.1.6): Suppose that we have the pattern (6, 6, 6). Then reordering indices 
according to the polygons in the pattern, we have by the estimates of Lemma 5.2 

(t) (6) (6) (6) (t-18) 
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Similarly, if we have the pattern (6, 6, 4), then we find 



^cTj < -1.52 pt- 1.52pt+0.S3pt+ {26 - 16) pt < 8 pt. 

(t) 



(Proof of 5.1.7): We use the same method as in the proof of 5.1.6. If we have the 
pattern (6, 6), and if t < 22, then 



Similarly, if we have the pattern (6,4, 4), then J2{t) ^^^^ than —1.52 pt+0. 33 pt+ 
0.33 pt+ 8pt< 8pt. 

(Proof of 5.1.8): Let the two adjacent degree six vertices be vi and V2- Let the 
third be ^3. The six triangles in the hexagon around give less than — 1.52 pt. 
The ten triangles in the hexagons around vi and f 2 give at most 



by the argument described in the case t = 26 of 5.1.2 (see Lemma 5.2). 

Suppose that t < 24. There remain at most 24 — 16 = 8 triangles, and they give a 
combined score of at most 8pt. The total score is then less than (—1.52 + 8) pt, as 
desired. 

Now assume that t = 26. Suppose there is a vertex v that is not adjacent to any 
of Vi, V2, or Vs- As in the proof of 5.1.8, the ten triangles in the two overlapping 
hexagons give less than pt. The other hexagon gives less than —1.52 pt. By Lemma 
5.2, the n triangles around v fall short of n points by at least (5 — 4.52) pt = 0.48 pt. 
Each of the remaining triangles gives at most Ipt. The score is then less than 
(0 - 1.52 + (26 - 16) - 0.48) pt = 8pt, as desired. 

(Proof of 5.1.9): This follows directly from the construction of the triangulation 
and the close-neighbor restrictions on the lengths of the edges of a quasi-regular 
tetrahedron. The lengths are between 0.8 < 2 arcsin(l/2.51) and 2 arcsin(2.51/4) < 
1.36. This completes the proof of the proposition. I 



ai < -1.52 pt- 1.52 pt+ {22 - 12) pt< 8pt. 



(t) (6) (6) (t-12) 




Section 6. Combinatorics 



Theorem 6.1. Suppose that a triangulation satisfies Proposition 4.2 and Proper- 
ties 1-9 of Proposition 5.1. Then it must the triangulation of Diagram 6.2 with 14 
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Diagram 6.2 













Proof: Fix a polygon centered at a vertex uq, such as the hexagon in Diagram 6.3. 
The six vertices Vi, . . . of the hexagon are distinct, for otherwise two distinct 
geodesic arcs on the sphere would run between uq and Vi for some i. This is 
impossible, because uq and Vi are not antipodal by Property 5.1.9. Similarly, the 
vertices of every other polygon of the triangulation are distinct. 



Diagram 6.3 




We then extend the polygon to a second layer of triangles. Each of the vertices Vi 
has degree four, five, or six, and two degree four vertices cannot be adjacent. The 
new vertices are denoted wi, . . . ,Wk- One example is shown in Diagram 6.4. 



Diagram 6.4 












As \/uo \'^2__^ - 




W5 ^- ^ -i 












W4 V A/^W3 
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There can be no identification of a vertex Vi with a vertex Wj, for otherwise there 
is a triangle (say with vertices Wj,Vk,uo,Vi = Wj) that is subtriangulated, con- 
trary to Property 5.1.3. Similarly, there is no identification of two vertices Wi and 
Wj, for otherwise it can be checked that there is a quadrilateral (with vertices 
Wi,Vk,uo,V£,Wj = Wi) that encloses more than one vertex, which is impossible 
by Proposition 4.2. A purely combinatorial problem remains. It is solved in the 
appendix. I 



Proposition 6.5. The triangulation of Theorem 6.1 scores less than 8pt. 

Proof: Our initial bound on the score comes by viewing the triangulation as made 
up of two hexagons and twelve additional triangles. By Lemma 5.2 and Calculation 
9.1, 

(6.6) - Xl^^ + X!^^ + -1.52 pt- 1.52 pt+ 12 pt^ 8.96 pt. 

(24) (6) (6) (12) 

A refinement is required in order to lower the upper bound to 8pt. 

If a vertex v has height l^;! > 2.2, then by Calculation 9.2, we find ai < 0.5 pt for the 
tetrahedra at v. Thus, in the Inequality 6.6, if the vertex v of some pentagon has 
height \v\ > 2.2, then the term ^(12) Pi ^^^J be replaced by ^(9) pt+ ^(3) 0-5 pt 
(there are three triangles in the pentagon that do not belong to the hexagon) , and 
the upper bound on the score falls to 7.46 pt. 

If a vertex of degree six has height \v\ > 2.05, then we claim that ^^^g^ < — 3.04 pt. 
In fact, by Calculation 9.7, the hexagon gives 

J^cTi < ^(0.389195 dih(^,) -0.435643) = 27r(0.389195)- 6(0.435643) < -S.04 pt. 

(6) (6) 

Estimate 6.6 is improved to 

^(Ti < -S.OApt- 1.52 pt + 12 pt= 7. Upt. 

(24) 

If the twelve tetrahedra have combined solid angle less than 6.48, then by Calcula- 
tion 9.9, 

^C7i < ^(0.446634 sol(,Si) - 0.190249) 

(12) (12) 

< 6.48(0.446634) + 12(-0.190249) < 11.039 

TVipn Flstimat.fi fi fi is imnrnvprl tn tViP bniinrl — 1 Ft'?, nf. — 1 Ft'?, nf.-^-^ 1 fl.'^Q nf. = 7 PQQ nf. 
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Now assume, on the other hand, that the combined sohd angle of the two hexagons is 
atmost47r-6.48. Set K = (47r-6.48)/12 and define ct' (5) := a{S) + {K-sol{S))/3. 
The sohd angle of one of the two hexagons is at most 6K. For that hexagon, we 
have 

^'{Si) = 5^ + (6K - 5^ so\{Si))/3 > ^i- 

(6) (6) (6) (6) 

By our previous estimates, we now assume without loss of generality that the heights 
l^;! of the vertices of triangles in the hexagon are at most 2.05, 2.2, and 2.2, the 
bound of 2.05 occurring at the center of the hexagon. By Calculation 9.15, 

Y^'i^i) < 0.564978 dihj -0.614725 = 27r(0.564978) + 6(-0.614725) < -2.5pt 

(6) (6) 

Estimate 6.6 becomes 

(24) (6) (6) (12) 

This completes the proof of the proposition. I 
Section 7. The Method of Subdivision 



The rest of this paper is devoted to the verification of the inequalities that have 
been used in Sections 5 and 6. In this section we describe the method used to 
obtain several of our bounds. We call it the method of subdivision. Let p{x) = 
Y^jCjx^ be a polynomial with real coefficients c/, where / = (ii, . . . ,in), x — 
(xi, . . . ,Xn) G M", and = x]^ ■ ■ •x'^n- It is clear that if C is the product of 
intervals [ai, 6i] x • • -x [a^, 6^] C in the positive orthant (oj > 0, for i = 1, . . . , n), 
then 

Mx e C, 

Pmin 

(C) < p{x) < (C), 

where 

]3mm(C) = X '^I^^ + X '^^^^ Pmax(C') = ^ ^/^^^ + ^ ^^^^ ' 

ci>0 c/<0 cj>0 cj<0 

Another bound comes from the Taylor polynomial p{x) = ^ (i/(a; — a)^ at a: 
(7.1) do+ X diib-aY <p{x) <do+ ^ diib-aY- 

di<0,1^0 di>0,IjtO 

If r{x) = p{x)/q{x) is a rational function, and if qmin{C) > 0, then 

w ^ PminiC) , / X ^ Pmax(C') 

e C, — — — r- < r{x) < 



qi{C,p) 92(c,p) ' 
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where qi{C,p) (resp. q2{C,p)) is defined as gmax(C') whenever Puiin{C) > (resp. 
Pma.x{C) < 0) and as qmin{C) otherwise. 

Let us define a cell to be a product of intervals in the positive orthant of M". By 
covering a region with a sufficiently fine collection of cells, various inequalities of 
rational functions are easily established. To prove an inequality of rational functions 
with positive denominators (say ri{x) < r2{x), for all x G C), we cover C with a 
finite number of cells and compare the upper bound of ri{x) with the lower bound 
of r2{x) on each cell. If it turns out that some of the cells give too coarse a bound, 
then we subdivide each of the delinquent cells into a number of smaller cells and 
repeat the process. If at some stage we succeed in covering the original region C 
with cells on which the upper bound of ri{x) is less than the lower bound of r2{x), 
the inequality is established. 

A refinement of this approach applies the method to the partial derivatives. If, for 
instance, we establish by the method of subdivision that for some i 

|^(x)>0, VxeC, 
dxi 

then we may compute an upper bound of p by applying the method of subdivision to 
the polynomial obtained from p by the specialization Xi = bi, where bi is the upper 
bound of Xi on C. Thus, we obtain an upper bound on a polynomial by fixing all the 
variables that are known to have partial derivatives of fixed sign and then applying 
the method subdivision to the resulting polynomial. Similar considerations apply 
to lower bounds and to rational functions with positive denominators. 

It is a fortunate circumstance that many of the polynomials we encounter in sphere 
packings are quadratic in each variable with negative leading coefficient (u, p, A, x 
in the next section). In this case, the lower bound is attained at a corner of the cell. 
Of course, the maximum of a quadratic function with negative leading coefficient 
is also elementary: —a{x — xq)^ + P < P, if a > 0. 

Section 8. Explicit Formulas for Compression, Volume, and Angle 

Many of the formulas in this section are classical. One can typically find them 
in nineteenth century primers on solid geometry. The formula for solid angles, 
for example, is due to Euler and Lagrange. For anyone equipped with symbolic 
algebra software, the verifications are elementary, so we omit many of the details. 
All formulas in this section will be valid for Delaunay simplices and for quasi-regular 
tetrahedra, unless otherwise noted. 

8.1. The Volume of a Simplex 
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As in the previous section, a cell in R" is a product of intervals in R". We define 
a function A : [4, 16]^ c R^ ^ R by 



Afxi , . . . , Xf\) = 


X^XA{—X^ -\- X') -\- X-\ — Xa. -\- X<^ -\- X(\\ 

+ a;2iC5(a^i - 3^2 + a^s + - a^s + xq) 
+ a;3a;6(a;i + 0:2 - 0:3 + 0:4 + ^5 - a^e) 


(8.1.1) 


— 0:23^33^4 — X\X->,X^ — X1X2XQ — X4X5XQ. 


We set Hi = s/xit for i = 1, . . 


. , 6. This relationship between Xi and i/i remains in 
inciex me etigeis oi a ounpiex aa ui uiagiani o.i.z. 




Diagram 8.1.2 




/\ ^\^4 




2 / \ 








\ 


We also define u : [4, 16]^ ^ 1 


^by 



u{xi,X2,xq) = {yi + y2 + VeMvi + 2/2 - y6)(yi - 2/2 + y6)(-yi + 2/2 + ye) 
(8.1.3) = —x^ — x^ — Xq + 2xiXq + 2xiX2 + 2x2Xq. 



Lemma 8.1.4. There exists a simplex of positive volume with edges of length 
Hi, ■ ■ ■ ,2/6 if cind only if A(a;i, . . . ,xq) > 0. If these conditions hold, then the 
simplex has volume A{xi, . . . , xq)^ /12. 

Proof: The function u = u{xi, X2,Xq) is quadratic in each variable, with negative 
leading coefficient, so the minimum of u, which is 0, is attained at a vertex of the 
cube [4, 16]"^. At the vertices where the minimum is attained, A(a;i, . . . ,xq) < 0. 

Assume A > 0. Then tt > 0, and a simplex exists with vertices 0, X = (j/i,0,0), 
Y = ^{*,u^ /yi,0), Z = {*,*,{A/u)^). Conversely, if the simplex exists, then it 
must be of the given form, up to an orthogonal transformation, so ti > (otherwise 
O, X, and Y are colinear) and A > 0. The volume is | det(X, Y, Z)\/6 = /12. 1 

Let C be a cell contained in [4, 16]^. The minimum of A on C is attained at a vertex 
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coefficient. To obtain an upper bound on A on a cell, we may use the method of 
Section 7. The restriction of A to the zero set of di^jdxx is 

4X4 

This is an upper bound on A(a;i, . . . ,xe). 
8.2. The Circumradius of a Simplex 

The circumradius of a face with edges y^, y^, and ye is 

v{y4,y5,y6) = 



u{yl,yl,yiy/^' 
Define p : [4, 16]^ ^ M by 



/ \ / \ oooooo 

+2X1X2X4X5+2X1XSX4XQ+2X2XSX5XQ. 



Lemma 8.2.2. Suppose A > 0; then p > 0. 

Proof: If A > 0, then the simplex of Lemma 8.1.4 exists. Let the coordinates 
of the circumcenter of the simplex be {x,y,z). Direct calculation shows that < 
x^ + y^ + z'' = p/{4A). I 

Corollary 8.2.3. Let 5 be a Delaunay simplex of positive volume with edges of 
lengths yi = ^/xi. The circumradius is |(p/A)2. 

Set 

x{xi,X2, X3, X4, X5, xe) = X1X4X5 + X1XQX4 + X2XGX5 + X2X4X5 + X5X3X6 

+ X3X4XQ — 2X5X^X4 — Xix\ — X2x\ — XsXq. 

We have 

p 2 . . . ,X6)2 

^-viy.,y.,ye) = 4^,^,,^,^^)^- 

The vanishing of x is the condition for the circumcenter of the simplex to lie in 
the plane through the face T bounded by the fourth, fifth, and sixth edges, x is 
positive if the circumcenter of S and the vertex of S opposite T lie on the same 
side of the plane through T and negative if they lie on opposite sides of the plane. 
We will say that T has positive orientation when x > 0- 



8.2.4. The function p is quadratic in each of the variables xi, 



, xq with negative 
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The derivative is A'^d{p/ A)/dxi — xi^b: ^6: ^i: ^2, xs, X4)x{xi, X2, xs, X4, x^, xq). 
This leads to bounds on the circumradius by the method of Section 7. 

The circumradius of a quasi-regular tetrahedron is increasing in xi, . . . ,xq if the 
orientation of each face is positive. When a face fails to have positive orientation, 
it satisfies the constraints of Section 3 (for example yi, 1/2, 1/3 G [2, 2.15], 1/4, 1/5, G 
[2.3,2.51], ri{y4,y5,yQ) > \/2, if it is the face opposite the origin). This allows us 
to determine bounds on the circumradius for most of the quasi-regular tetrahedra 
we encounter in this paper by inspection. For example, in Section 9 we will study 
the simplices constrained by yi G [2,2.1]. An upper bound on the circumradius is 
rad(,S(2.1,2.1,2.1,2.1,2.1,2.1)). 

8.2.5. Proof of Lemma 3.4. In the notation of Section 3.4, let T be the given face, 
and let Vq be a vertex satisfying the conditions of the lemma. Let S be the simplex 
at the origin vq — 0, whose first, second, and third edges abut a-t Vi, for i — 1,2,3. 
Suppose for a contradiction that \vo — vs\'^ = xs > 2.15^. In light of the results of 
this section, a contradiction follows if xi^i, ■ ■ ■ ,X6) > 0, for Xi > 4 for z = 1,2, 
X3 > 2.15^, 2.3^ <Xi< 2.51^ for i = 4,5,6. (The lower bound of 2.3 comes from 
Remarks 3.2 and 3.3: the circumradius of T is at least -\/2). Since T is acute. 




= X4{-X4 -1-0:5-1- xe) > 0. 



Similarly, x is increasing in X2 and x^. Thus, to minimize x, we take xi = X2 = 4, 
xs = 2.15^. But X is quadratic in each variable X4, x^, and xq with negative 
leading coefficient, so the minimum occurs at a corner point of [2.3^,2.51^]"^. We 
find X > X(2^ 22, 2.152, 2.512, 2.51^, 2.51^) 0.885 > 0. I 



For the set of Delaunay stars to be a compact topological space, simplices of zero 
volume must be included. Since the circumradius remains bounded, the degen- 
erate Delaunay simplices of zero volume are planar quadrilaterals that possess a 
circumscribing circle. 

8.3. Dihedral Angles 

Let dih(5') be the dihedral angle of a simplex S along the first edge. It is the 
(interior) angle formed by faces with edges (2, 1, 6) and (1, 3, 5). 

Lemma 8.3.1. 

dA/dx4 



cosdih(S') 



{u{xi, X2, xe)u{xi, X3, X5)y/'^ ' 
where u is the function defined by Equation 8.1.3. 
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The partial derivative of cosdih(S') with respect to is 

2a;i(aA/ax2) 

u{xi,X2, X6y/'^u{xi,X3, 0:5)3/2 ' 

SO that the sign of the partial derivative is determined by the sign of OA/ 8x2- 
Similar considerations apply to the partial derivatives with respect to X2, X5, and 
xe- 

Lemma 8.3.2. Let 5" be a quasi-regular tetrahedron. Then the dihedral angle of 
S along the any edge is at most arccos ( "95999^ ) — dihmax ~ 1.874444. 

Proof: Suppose that dih(S') > 7r/2, so that the numerator n = OA/ 8x4 in Lemma 
8.3.1 is negative. To bound cosdih(S') from below, we minimize u{xi,X2,xq), 
u{xi, x^, X5), and the numerator n over the variables X2,X3,X5, and xq. 

We have diU = xj + Xk — Xi > on the indicated domain, for i = 2, 3, 5, and 6, so 
that 

u{xi, X2, xq)u{xi, xs, 0:5) > u{xi, 2^, 2^)^ = {-x\ + 16x1)^. 
Similarly din = xi-\- xj — x^ > 0, for i = 2, 3, 5, and 6, so that 

n(xi, . . . , Xq) > n(xi, 2^, 2^, X4, 2^, 2^) = xi(16 — xi — 2x4). 

Thus, 

^•^,^c^ ^ 16-0:1-20:4 ^ 20:4 2(2.51)2 -29003 

> cosdih(5) > = 1 > 1 ^ — . 

^ ^- 16-0:1 16-0:1- 16- 2.512 96999 



8.4. The Solid Angles of a Delaunay Simplex 



Let S* be a Delaunay simplex with vertices vq, ■ ■ ■ , t's- By the solid angle soli(S') of 
the simplex at the vertex i, we mean 3vol(5' flSj), where Bi is a unit ball centered at 
Vi. For simplicity, suppose that vq is located at the origin, and that in the notation 
of Diagram 8.1.2, the vertices vi,V2, are the edges of lengths yi, y2, and y^. We 
write sol for soIq. Set 

(8.4.1) ^ ^ ^ 

a{yi, 2/2, • • • , ye) = ymys+^yi {yl +yl-yl) + ^y2{yl+yl -yl) + -^y^iyl+yl-yl)- 



Lemma 8.4.2. 

2a 

sol(5') = 2 arccot(^Yy^). 
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Proof: (See [H2, p.64]). We use the branch of arccot taking values in [0,7r]. I 

The function a is increasing in yi, 1/2, and ys on [2,4] and is decreasing in the 
variables 7/4, j/s, and uq on the same interval. 

We claim that a > for any Delaunay simplex. Wc will prove this in the case 
A > 0, leaving the degenerate case A = to the reader. The area of the face T 
opposite the origin is at most 4\/3 (since its edges are all at most 4). There exists 
a plane that is tangent to the unit sphere between this face and the unit sphere 
[HI, 2.1]. The area of the radial projection of T to the this plane is less than the 
area of a disk D of radius 1.5 on the plane centered at the point of tangency. The 
solid angle (that is, the area of the radial projection of T to the unit sphere) is 
less than the area of the radial projection of D to the unit sphere. This area is 
27r(l — cos(arctan(1.5))) < tt. Lemma 8.4.2 now gives the result. This allows us to 
use Lemma 8.4.2 in the form sol(5') = 2 arctan(A-'^/^/(2a)). 

8.4.3. The solid angle is the area of a spherical triangle. Let x, j/, and z be the 
cosines of the radian lengths of the edges of the triangle. By the spherical law of 
cosines, the solid angle, expressed as a function of a;, j/, and z is 



The partial derivative with respect to x of this expression for the solid angle is 
{—l — x + y + z)/{{x + l)\/t)^ where t = l — x^ — y'^ — z^ + 2xyz. The second derivative 
of this expression evaluated at —1 — x + y + z = is —2(1 — y){l — z)t~^/'^ < 0. So 
the unique critical point is always a local maximum. If neither edge at a vertex of 
a spherical triangle is constrained, then the triangle can be contracted or expanded 
by moving the vertex. If the lengths of the edges are constrained to lie in a product 
of intervals, then the minimum area occurs when two of the edges are as short as 
possible and the third is at one of the extremes. The maximum area is attained 
when two of the edges are as long as possible and the third is at one of the extremes 
or at a critical point: x = y + z — 1 (or it symmetries, y = x + z — 1, z = x + y — 1). 

8.5. The Compression of a Delaunay Simplex 

The compression T{S) of a Delaunay simplex S is defined as 



c{x, y, z) + c{y, z, x) + c{z, x, y) — tt, c(a;, y, z) := arccos( 



X — yz 



V(i-y2)(i-^2) 



4 



(8.5.1) 




where 



(8.5.2) 



37r + 12arccos(l/V3) 



0.72. 
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We let ao = (1(1/1,1/2,1/3,1/4:11/511/6): where a is the function defined in 8.4.1. We 
also let ai, a2, and be the functions a for the vertices denoted vi, V2, and V3 in 
Diagram 8.1.2. For example, ai = j/s, J/B; J/4) 2/2) J/s)- Set t = y/A/2. By the 
results of Sections 8.1 and 8.4, we find that r(5') = —5oct i/^+Yloi'^/^) arctan(t/ai). 
Recall that > on a cell C. We wish to give an elementary upper bound of F 
on C. Set, for t > and a = (oq, ai, a2, 03) e M^, 

^ t 2 

7(^) a) := — ^ + o arctan( — ). 

The partial derivative of 7 with respect to is — 2t/(3(t^ + af)) < 0, so 7(t, a) < 
7(t, a~), where a~ = (a^, aj~, a^, 03 ) is a lower bound of a on C, as determined by 
Section 8.4. We study ^(t) = ^{t, a~) as a function of t to obtain an upper bound 
on V{S). We note that 

VW = 4^ + 1 1 ^ and y'(0 = -I E ^ 0- 

Upper and lower bounds on t are known from Section 8.1. An upper bound on 7 
for t in [tmiiD^max] IS the maximum of -^(tmin) and -^(tmax)) where £ is the tangent 
to 7 at any point in [tmiiD ^max]- 



8.6. Voronoi cells 

We assume in this section that the simplex S has the property that the circumcenter 
of each of the three faces with vertex at the origin lies in the cone at the origin 
over the face. This condition is automatically satisfied if these three faces are acute 
triangles. In particular, it is satisfied for a quasi-regular tetrahedron. 

When the cone over a quasi-regular tetrahedron S contains the circumcenter of S, 
Section 2 sets vor(S') = — 4(5octvol(S'o) -|- 4sol(5')/3, where is the intersection of 
5' with a Voronoi cell at the origin. Otherwise, vor(S') is defined as an analytic 
continuation. This section gives formulas for vol(S'o). 

As usual, set S = S{yi, . . . , j/e) and Xi = yf. Suppose at first that the circumcenter 
of S is contained in the cone over S. The polyhedron Sq breaks into six pieces, 
called the Rogers simplices. A Rogers simplex is the convex hull of the origin, the 
midpoint of an edge (the first, second or third edge), the circumcenter of a face 
along the given edge, and the circumcenter of S. Each Rogers simplex has the form 

R. = B.(a. h. r^ := S(o,. h. c. ic^ - h^)^/\ ic^ - n^)^/''. (h^ - a^)^/^) 
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for some 1 < a < b < c. Here a is the half-length of an edge, b is the circumradius 
of a face, and c is the circumradius of the Delaunay simplex. 

The volume is vol(i?(a, &, c)) = a{b'^ - 0^)1/2(^2 _ b^y/^/6. The density 5{a, b, c) of 
R — R{a, 6, c) is defined as the ratio of the volume of the intersection of R with a 
unit ball at the origin to the volume of R. It follows from the definitions that 

(8.6.1) vot{S) = J]4vol(i?(a,6,c))(-5oct + 5(a,6,c)), 

where c is the circumradius rad(S'), and (a, b) runs over the six pairs 

(yi/2, viVi, 2/2, ye)), (2/2/2, r]{yi, y2, ye)), 
(^2/2, v{y2, ys, 2/4)), (2/3/2, r]{y2, 2/3, 2/4)), 
(2/3/2, 77(2/3, 2/1,2/5)), (2/1/2, 7/(2/3, 2/1, 2/5))- 

Upper and lower bounds on vol( J?(a, 6, c)) follow without difficulty from upper and 
lower bounds on a, 6, and c. Thus, an upper bound on 5(a, 6, c) leads to an upper 
bound on vor (5). The next lemma, which is due to Rogers, gives a good upper 
bound [R]. 

Lemma 8.6.2. The density 5{a, b, c) is monotonically decreasing in each variable 
for 1 < a < 6 < c. 

Proof: Let 1 < ai < 61 < ci, 1 < 02 < 62 < C2, 02 < ai, 62 < ^1, C2 < ci. The 
points of R{ai, 61, Ci) arc realized geometrically by linear combinations 

si = Ai(ai, 0, 0) + A2(ai, (6? - a?)^/^, 0) + A3(ai, (6? - al)'/', (cf - 62)i/2)^ 

where Ai, A2, A3 > and Ai + A2 + A3 < 1. The points of i?(a2, 62, C2) are realized 
geometrically by linear combinations 

S2 = Ai(a2, 0, 0) + A2(a2, {bj - 4)'/^ 0) + A3(a2, {bl - alfl\ [cl - blf'^), 

with the same restrictions on Aj. Then 

|si|2 - |S2|2 = Ai(Ai + 2A2 + 2A3)(a2 - al) + A2(A2 + 2\-i){b\ - b\) + \l{c\ - cl). 

So |sip > |s2p. This means that the linear transformation Si 1— > S2 that carries the 
simplex to 5*2 moves points of the simplex 5'i closer to the origin. In particular, 
the linear transformation carries the part in 5*1 of the unit ball at the origin into the 
unit ball. This means that the density of R{a\, bi,ci) is at most that of R{a2, &2, C2)- 



If the circumcenter of S is not in the cone over S, then the analytic continuation 
gives 

vor(S') = ^4:eRVo\{R{a,b,c)){-5oct + S{a,b, c)). 
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where cr — 1 if the face of the Delaunay simplex S corresponding to R has positive 
orientation, and eji = —1 otherwise. (The face of S "corresponding" to R{a, h, c) is 
the one used to compute the circumradius h.) 



8.6.3. A calculation based on the explicit coordinates of S and its circumcenter 
given in 8.1.4 shows that 

eHVol(i?(,,/2,ry(,„,„,e),rad(5))) = 48.(.„ .e)A(.„ . . . , .e) ' 

If the circumcenter of 5' is not contained in 5', then the same formula holds by- 
analytic continuation. By definition, en = —1 exactly when the function x is 
negative. Although this formula is more explicit than the earlier formula, it tends 
to give weaker estimates of vor(S') and was not used in the calculations in Section 
9. 



8.6.4. There is another approximation to vor(5') that will be useful. Set Sy = 
S'(2, 2, 2, J/, y, J/). (We hope there is no confusion with the previous notation Si.) 
For 1 < a < 6 < c, let vol{R{a, b, c)) = a{{b'^ - a'^){c^ - b'^))^/^/6 be as above. Set 
r(a) = vol(i?(a,r/(2,2,2a),1.41)). 

Lemma 8.6.5. Assume that the circumradius of a quasi-regular tetrahedron S is 
at least 1.41, and that 6 < yi + 1/2 + Us ^ 6.3. Set a — {yi + 1/2 + Us — 4)/2. Pick y 
to satisfy sol{Sy) = sol(S'). Then 

vot{S) < voT{Sy) - 85oct{l - l/a^ia) 



Proof: If S is any quasi-regular tetrahedron, let S'tan be the simplex defining 
the "tangent" Voronoi cell, that is, S'tan is the simplex with the same origin that 
cuts out the same spherical triangle as S on the unit sphere, but that satisfies 
j/i = j/2 = Z/3 = 2. The lengths of the fourth, fifth, and sixth edges of S'tan are 
between ^/S — 2(2.3) = VSA and 2.51. The faces of ^tan are acute triangles. A 
calculation similar to the proof in 8.2.5, based on x(3-4, 3.4, 4, 4, 4, 2.51^) > 0, shows 
that the circumcenter of S'tan is contained in the cone over Stan- 

Since 5'tan is obtained by "truncating" S, we observe that vor(5') = vor(S'tan) — 
45octvol(5'\5'tan), whcrc S and S'tan are the pieces of Voronoi cells denoted Sq in 
Section 2 for S and 5'tan, respectively. By a convexity result of L. Fejes Toth, 
vor(S'tan) < vor(S'y) [FT,p.l25]. (This inequality relies on the fact that the circum- 
center of 5'tan is contained in the cone over -Stan-) 

Let ai be the half-length of the first, second, or third edge. Since vol(-R(ai, 6, c)) is 

inrrfiflsincr in r wfi nht.a.in a. Inwfir hniinri nn \rn]( R( n.-i h r\\ fnr = 1 41 tlip Inwpr 
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bound on the circumradius of S. The function vol(i?(ai, b, 1-41)), considered as a 
function of 6, has at most one critical point in [ai, 1.41] and it is always a local 
maximum (by a second derivative test). Thus, 

(8.6.6) vol(i?(ai, b, c)) > min(vol(i?(ai, 6„,i„, 1.41)), vol(i?(ai, b^^^, 1.41))), 

where 6min and 6max are upper and lower bounds on 6 e [ai, 1.41]. A lower bound 
on b is 77(2, 2, 2ai), which means that vol(i2(ai, femiiu 1-41)) may be replaced with 
r(ai) in the Inequality 8.6.6. By Heron's formula, the function rj{a, b, c), for acute 
triangles, is convex in pairs of variables: 



VaaVbb - vlb = '' ' ^HuL^^ ' " ' " ' " ' " ^ > 0. 



7j%a^ + 62 - c2)(a2 - 62 + c^){-a ^ + b^ + c^){a^ + 6^ + c') 



Thus, an upper bound on 6max is ^7(2, 2.51, 2a), where a = (j/i + 2/2 + J/3 — 4)/2. 
This means that 

vol(i?((ai,6o,ax,1.41))) 

may be replaced with the function vol(i?(ai, 77(2, 2.51, 2a), 1.41)) in the Inequality 
8.6.6. Now 1 < ai < a, and vol(i?(ai, 6, c)) is decreasing in the first variable (for 
1 < ai and 6 < ^), so we may use the lower bound 

vol(i?(a, 7/(2, 2.51, 2a), 1.41)) 

instead. By Calculations 9.20.2 and 9.20.3, wc conclude that vol(f2(ai, 6, c)) > r(a). 
This lower bound is valid for each Rogers simplex. The volume of S'\S'tan is then 
at least 

((1 - 8/y?) + (1 - 8/y|) + (1 - 8/y|))2r(a). 
The concavity of 1 — 8/y"^ gives 

Ei=i(l - ^/vl) > 1 - We have established 

that 

vol(^\^tan) > (1 - l/a3)2r(a). 
The result follows. I 



8.6.7. We conclude our discussion of Voronoi cells with a few additional comments 
about the case in which analytic continuation is used to define vor(5'), with 5" 
a quasi-regular tetrahedron. Assume the circumcenter c of lies outside S and 
that the face T of S with negative orientation is the one bounded by the first, 
second, and sixth edges. It follows from Section 3 that 7/1,2/2,2/6 ^ [2.3,2.51] and 
2/3,2/4,2/5 e [2,2.15]. 

Let pi (resp. P2) be the point on T equidistant from the origin, v^, and vi (resp. 

|2 _ ^1 Xiu{xi,X2jX6){-Xi + Xs + X^)'^ 
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ei? = — 1 for the Rogers simplex with vertices the origin, po, c and vi/2. It hes 
outside S. The other Rogers simplex along the first edge has ej? = 1, so that the 
part common to both of these Rogers simplices cancels in the definition of vor(S'). 
This means that vor(5') becomes the sum of the usual contributions from the two 
Rogers simplices along the third edge, 

4vol(i?(a, 6, c))(-5oct + ^(a, c)) 

for (a,6,c) = (yi/2, ?/3, ys), bi|)) and (2/2/2, 77(2/2, 2/3, 2/4), |P2|), and 

v{^S") := 45oetVol(^") - 4sol(5")/3, 

where S" is the convex hull of 0, po, Pi, P2, and c. 

We claim that in any cell satisfying the constraints given above, \p\\ is minimized 
by making yi, 2/2, 2/6 as large as possible and 2/3, 2/4, 2/5 as small as possible. To 
show this, one has to write out the derivatives explicitly from the formula for 
given above. The partial derivatives with respect to x^-, ■ ■ ■ factor into 
products of the polynomials d/!:^/dxi, for z = 3,4, 5, a^i, {x:i + x^ — a^i), it(xi, X2, a^e), 
{—x\X4^-\-X2X^ + x\Xq — x^xq)., and {xiX2—X2X:i—xiX4^-\-x:iXQ). The signs of these 
polynomials are easily determined. The partial with respect to xi is complicated 
(the numerator has 88 terms), and we had to resort to the method of subdivision 
to determine its sign. We omit the details. 

To complete our estimate, we describe an upper bound on v{S"). The ratio 
sol(^")/(3vol(5")) is at least l/rad(^)3, because S" is contained in a sphere of 
radius rad(/S'), centered at the origin. We have 

vol(S") 1, 

— ^ — < g|Po -Pil \c-po\ \po\, 
because the convex hull of 0,po,Pi, and c, which is one side of S" , is a pyramid 
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course, |c — PoP = rad(S')^ — ?7(|/i, I/2, ye)^- We now have a bound on v{S") in 
terms of quantities that have been studied in Section 8, if we rely on the bound 
bo — Pi\ < 0.1381. Write pi = Pi{yi, ■ ■ ■ , j/e)- This bound is obtained from the 
following inequalities. 

bo -Pi\ = ivivi, y2, yef - - {\pif - 

< {riiyi, 2.51, 2.51)^ - a;i/4)i/2 _ (|^^(^^^ 2.51, 2, 2, 2, 2.51)1^ - xi/^y/^ 

< (77(2.51, 2.51, 2.51)^ - 2.5lV4)^/^- 

(bi(2.51, 2.51, 2, 2, 2, 2.51)p - 2.5lV4)^/^ 

< 0.1381. 

The inequality that replaces yi with 2.51 results from Calculation 9.21. 
8.7. A Final Reduction 

Let S he a, Delaunay simplex. Suppose that the lengths of the edges yi, j/s, and ye 
are greater than 2. Let S' be a simplex formed by contracting the vertex joining 
edges 1, 5, and 6 along the first edge by a small amount. We assume that the 
lengths y[, 1/5, and y'^ of the new edges are still at least 2 and that the circumradius 
of S' is at most 2, so that S' is a Delaunay simplex. 

Proposition 8.7.1. T{S') > T{S). 

We write sol^, for i = 1, 2, 3, for the solid angles at the three vertices pi, p2, and ps 
of S terminating the edges 1,2, and 3. Let p'l be the vertex terminating edge 1 of 
S'. Similarly, we write sol^, for z = 1, 2, 3, for the solid angles at the corresponding 
vertices of S'. We set vol(F) = vol(5') — vol{S') and Wi = solj — sol^. It follows 
directly from the construction of 5" that W2 and are positive. The dihedral 
angle a along the first edge is the same for S and S'. The angle (3i of the triangle 
(0,pi,Pi) at pi is less than the angle of the triangle {0,p[,pi) at p'l, for z = 2, 3. 
It follows that wi is negative, since —wi is the area of the quadrilateral region of 
Diagram 8.7.2 on the unit sphere. 
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Lemma 8.7.3. Soct^oliV) > t(;2/3 + tus/S. 

The lemma immediately implies the proposition because wi < and 
V{S') - V{S) = -wi/Z - W2/?> - «;3/3 + 5octvol(F). 



Proof: Let T' be the face of S' with vertices p'^, ^>2, and p^. We consider S" as 
a function of t, where t is the distance from pi to the plane containing T' . (See 
Diagram 8.7.2.) It is enough to establish the lemma for t infinitesimal. 

As shown in Diagram 8.7.2, let Vi be the pyramid formed by intersecting V with 
the plane through pi that meets the fifth edge at distance to = 1-15 from p3 and 
the sixth edge at distance to from p2. Also, let the intersection of V with a ball 
of radius to centered at pi be denoted V^, for i = 2, 3. For t sufficiently small, the 
region Vi is (essentially) disjoint from V2 and V3. 

We claim that vol(Vi) > vol(V2nV3). Let 9 be the angle of T' subtended by the fifth 
and sixth edges of S' . Then vol(Vi) = Bt/Z, where B is the area of the intersection 
of Vi and T': 

(8.7.4) S = sin^(y^-to)(2/6-^o)/2. 

Since S' is a Delaunay simplex, the estimates 7r/6 < 9 < 27r/3 from [HI, 2. 3] hold. 
In particular, sin^ > 0.5, so 5 > 0.25(2 - 1.15)^ and vol(Fi) > 0.06 1. 

If vol(V2 n V3) is nonempty, the fourth edge of 5" must have length less than 2to- 
The dihedral angle a' of V along the fourth edge is then less than the tangent of 
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the angle, which is at most t + in Landau's notation. As in [HI, 5], we obtain 

the estimate 

vol(F2 n Vs) < a' r tl -t^dt= fa-l)^(2to + lK ^ 0.025 1 + Oit^). 
Jl 3 

This estabhshes the claim. 
Thus, for t sufficiently small 

Soctyol{V) > Soct{yol{Vi) + Y0\{V2) + Y0l{Vs) - Vo\{V2 H F3)) 

>(^""o)(f.f)>1.09(f.f). 



Section 9. Floating-Point Calculations 

This section describes various inequalities that have been established by the method 
of subdivision on SUN workstations. The full source code (in C++) for these calcu- 
lations is available [H6]. 

Floating-point operations on computers are subject to round-off errors, making 
many machine computations unreliable. Methods of interval arithmetic give users 
control over round-off errors [Int] . These methods may be reliably implemented on 
machines that allow arithmetic with directed rounding, for example those conform- 
ing to the IEEE/ ANSI standard 754 [W],[IEEE], [NR]. 

Interval arithmetic produces an interval in the real line that is guaranteed to contain 
the result of an arithmetic operation. As the round-off errors accumulate, the 
interval grows wider, and the correct answer remains trapped in the interval. Apart 
from the risk of compiler errors and defective hardware, a bound established by 
interval arithmetic is as reliable as a result established by integer arithmetic on 
a computer. We have used interval arithmetic wherever computer precision is a 
potential issue (Calculations 9.1 - 9.19, in particular). 

Every inequality of this section has been reduced to a finite number of inequalities 
of the form r{xo) < 0, where r is a rational function of x e M"^ and xq is a given 
element in W^. To evaluate each rational expression, interval arithmetic is used to 
obtain an interval Y containing r{xo). The stronger inequality, y < for all y 
which may be verified by computer, implies that r{xo) < 0. 

To reduce the calculations to rational expressions r(xo), rational approximations to 
the functions ^/x, arctan(a:), and arccos(x) with explicit error bounds are required. 
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as TT, \/2, and 5oct) with explicit error bounds are also required. These were obtained 
in Mathematica and were double checked against Maple. 

Let S = S{yi,... jI/q) be a quasi-regular tetrahedron. We label the indices as 
in Diagram 8.1.2. Let F = r{S{yi, . . . ,2/6)) be the compression. We also let 
dih = dih{S{yi, ... ^ye)) be the dihedral angle along the first edge, and let sol = 
sol(S'(2/i, . . . , j/e)) be the solid angle at the origin, that is, the solid angle formed 
by the first, second, and third edges (j/i, j/25 J/a) of 5. 

All of the following inequalities are to be considered as inequalities of analytic func- 
tions of yi, . . . ,2/6. Although each of the calculations is expressed as an inequality 
between functions of six variables, Lemma 8.7.1 has been invoked repeatedly to 
reduce the number of variables to three or four. For instance, suppose that we wish 
to establish /(sol, F) < 0, where /(sol, F) is an expression in T{S{yi, ... ,2/6)) and 
sol(S'(yi, . . . ,ye)). Invoking 8.7.1 three times, we may assume that the vertices 
marked vi, V2, and vs in Diagram 8.1.2 each terminate an edge of minimal length. 
It is then sufficient to establish the inequality in seven situations of smaller dimen- 
sion; that is, we may assume the edges i e / have minimal length, where / is one 



Similarly, for an inequality in F and dihedral angle, we reduce to the seven cases 



The first two calculations are inequalities of the compression F of a quasi-regular 
tetrahedron. 

Calculation 9.1. T <lpt. 

This first inequality and Calculation 9.3 are the only ones that are not strict inequal- 
ities. Set So — S{2, 2, 2, 2, 2, 2). By definition, F(S'o) = 1 pt. We must give a direct 
proof that gives the maximum in an explicit neighborhood. Then we use the 
method of subdivision to bound F away from 1 pt outside the given neighborhood. 
An infinitesimal version of the following result is proved in [HI] . 



Lemma 9.1.1. If y, e [2, 2.06], for i = 1, . . . , 6, then r{S{yi, ... , ye)) < F(5o), 
with equality if and only if S{yi, . . . ,yQ) — Sq. 

Set aoo = a(2, 2, 2, 2, 2, 2) = 20, Aq = A{2^, 2\ 2^, 2^, 2\ 2^), = v^/2 = A^^2, 
and 60 = (2/3)(l + ^0/^00)"^ = 50/81. Set / = max,(2/, - 2) < 0.06 and a" = 



„/ooooi -Poi -Poi f\ on 1 o -P o -P2 A „ ;„ c — .j-; — o — j- j- . /a" /o 



of 



{1,4}, {2,5}, {3,6}, {4,5}, {4,6}, {5,6}, {1,2,3}. 



/ = {1,4}, {2,5}, {3,6}, {1,2,3}, {1,5,6}, {3,4,5}, {2,4,6}. 
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We have arctan(a;) < arctan(a;o) + {x — xo)/{l + Xq), if x^xq > 0. We will verify- 
below that A > Ao under the restrictions given above. Then 



3 



6 ^ \ai 000/ 
SoctjA - Aq) ^ t - to 

= ^(*°)- 24(t + t„) — 

+ -2- 2^ (aoo + 



an ^ ai 



i=0 " i=0 



< T{So) + ——— — — — + — + > (oqo - Oi) 

t + to \ 24 a-J alof^^^ 

3 



(9.1.2) +^J^(a„„-a.)^ 



i=0 



Set Co = -5oct/24 + 6o/a- > 0. Write = 2 + fi, with < fi < f. Set 
a;i = 4 + = = 4 + 4/, + /f. Then < 6^/4, for z = 1, . . . , 6. Set e = 4/ + J^. 
Set A(ei, . . . jCg) = A(xi, . . . jXe). We find that d'^A/de2dei = 4 + 64+65-66 > 
and similarly that A/ dcidei > 0, for i = 3,5,6. So dA/dei is minimized by 
taking 62 = 63 = 65 = ee = 0, and this partial derivative is at least 

dA , , y ^ 

— (ei, 0, 0, 64, 0, 0) = 16 - 8ei - 2eie4 - > 16 - 8e - 3e^ > 0. 
aei 

Thus A > Ao. The partial derivative dA/dei is at most 

dA . 22 

— — (ei, e, e, 64, e, e) = 16 + 16e — 8ei + 4ee4 — 2eie4 — 64 < 16 + 16e + 4ee4 — 64 
oei 

< 16 + 16e + 3e^ 
So A - Ao < (16 + 16e + 3e2)(ei + h eg). 

We expand X]i=o(oi — ooo) = 5(6i + ■ ■ - + 66) + /ii + /i2 + ^3 as a sum of homogeneous 
polynomials hi of degree i in /i, . . . , /e. A calculation shows that hi = 0. Also /i2 
is quadratic in each variable fi with negative leading coefficient, so /12 attains its 
minimum at an extreme point of the cube [0,/]^. A calculation then shows that 
h2 > —14/^ on [0, f]^. By discarding all the positive terms of hs, we find that 



hs > -flh - hfl - fih - /2/I - /I/6 - /3/I > -6/^ 
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Thus, ^j=o('^* ~ '^oo) > 5(ei + h ee) — 14/^ — 6/^. Since e = max(ei), we find 

that / < e/4 < (ei + h e6)/4. This gives 

^(aoo - ai) < (ei H h ee){-5 + ^ + -^)- 

i=0 



Similarly, we expand ^j=o('^oo ~ '^«)^ ^ ^ polynomial in /i, . . . , /e. To obtain an 
upper bound, we discard all the negative terms of the polynomial and evaluate all 
the positive terms at (/i, . . . , /e) = (/,.. . , /). This gives 

3 

XI ("00 -aif < 4944/2 + 7296/^ + 3684/^ + 828/^ + 73/^ 

i=0 

< ^(ei + • • • + e6)(4944/ + 7296/^ + 3684/^ + 828/^ + 73/^). 
We now insert these estimates back into Inequality 9.1.2. This gives 



f^^^l ~ 1^°^ < 1^(16 + 16e + ae^) + ^(-5 + 3.5/ + 1.5 f) 
[ei + \- ee) zto ago 



(9.1.3) + y°_ (4944/ + 7296/2 + 3684/^ + 828/^ + 73/^). 

The right-hand side of this inequality is a rational function of /. (Both a~ and e 
depend on /.) Each of the three terms on the right-hand side is increasing in /. 
Therefore, the right-hand side reaches its maximum at / = 0.06. Direct evaluation 
at / = 0.06 gives T{S) - T{So) < -0.00156(ei + • • • + eg). I 



To verify Calculation 9.14, the computer examined only 7 cells. Calculation 9.1 
required 1,899 cells. Calculation 9.6.1 required over 2 million cells. The number 
of cells required in the verification of the other inequalities falls between these 
extremes. 



Calculation 9.2. T < 0.5 pt, if yi e [2.2,2.51]. 



The next several calculations are concerned with the relationship between the di- 
hedral angle and the compression. 

Calculation 9.3. 

dih{S{yi, ... , ye)) > dih^in := dih(5(2, 2.51, 2, 2, 2.51, 2)) ^ 0.8639. 



Since this bound is realized by a simplex, we must carry out the appropriate local 

1 ;„ ;„ „ — 1 ] „f C/'o o ci o o o ci o^ 
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Lemma 9.3.1. Suppose that 2 < i/i < 2.2, for i = 1,3, 6, and that 2 < i/i < 2.51, 
for 1 = 2,4, 5. Then dih{S{yi, ... , ye)) > dih^in- 

Proof: This is an apphcation of Lemma 8.3.1. By that lemma, dih(S') is increasing 
in X4, so we fix X4 = 4. The sign of 9cosdih /dx2 is the sign of dA/dx^, and a 
simple estimate based on the explicit formulas of Section 8 shows that dA / 8x3 > 
under the given constraints. Thus, we minimize dih(S') by setting X2 = 2.51^. By 
symmetry, we set = 2.51^. 

Now consider dih(5') as a function of xi, X3, and xq. The sign of 9cosdih /dx^ is 
the sign of dA{xi, 2.51^, X3, 2^, 2.51^, XQ)/dx2- The maximum of this partial (about 
—2.39593) is attained when xi, x^, and are as large as possible: Xi = x^ = Xq = 
2.2^ . So dih(5') is increasing in X3. We take X3 = 4, and by symmetry xq = 4. 

We have 

acosdih(S') _ 2ti(xi) 

dxi u{xi,X2, X6)^/'^u{xi, X3, Xs)"^/^ ' 

where ti{xi) ^ 464.622 - 1865. 14xi + 326.954xf - 92.9817xf + 4xf. An estimate 
of the derivative of ti{x) shows that ti{x) attains its maximum at xi = 4, and 
ti(4) < —6691 < 0. Thus dih(5') is minimized when xi = 4. I 



Calculation 9.4. T < 0.378979 dih -0.410894. 

Lemma 9.5. If 5*1, S2, Ss, and 5*4 are any four tetrahedra such that dih(5'i) + 
dih(52) + dih(53) + dih(54) > 2n, then T{Si) + ■■■ + T{S4) < 0.33 pt. 

This lemma is a consequence of the following three calculations, each established 
by interval arithmetic and the method of subdivision. 

Calculation 9.5.1. 

r < -0.19145 dih +0.2910494, 
provided dih e [dih(5(2, 2, 2, 2, 2, 2)), 1.42068]. 

Calculation 9.5.2. 

r < -0.0965385 dih +0.1562106, 
provided dih e [1.42068, dih(S'(2, 2, 2, 2.51, 2, 2))]. 

Calculation 9.5.3. 

r < -0.19145 dih +0.31004, 

nrnvidpd dih > r[\h(S(9.^ 9.^ 9.^ 2 51 . 9^ 
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To deduce Lemma 9.5, we consider the piecewise linear bound ^ on F obtained from 
these estimates. The hnear pieces are £i{x) < 1 pt, £2, £3, and £4 on [dihmin, rfi], 
[di,d2], [d2,d3], and [d3,dih^^^] , where di = dih(5(2, 2, 2, 2, 2, 2)), ^2 = 1.42068, 
and ds = dih(5(2, 2, 2, 2.51, 2, 2)). (See Calculation 9.1, and Lemma 8.3.2.) Dia- 
gram 9.5.4 illustrates these linear bounds. (There are small discontinuities at di, 
d2, and ds that may be eliminated by replacing £i by £i + Ci, for some > 0, for 
i = 1,2, and 4.) We then ask for the maximum of i{ti) + £{t2) + ^(ts) + ^(^4) 
under the constraint {ti +^2+^3 +^4)/4 > 7^/2 G [(i2, ds]. Since £{x) is constant on 
(dihmin, (ii) and decreasing on [di , dihmax] , we may assume that U > di, for all i. 
Since the slope of £2 is equal to the slope of £4, we may assume that U > d2, for all i, 
or that ti < ds, for all i. Uti < ds, for all z, then we find that ti > 27r — Scia > (^2- So 
in either case, ti > d2, for all i. By convexity, an upper bound is 4:£s{7r/2) < 0.33 pt. 
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Lemma 9.6. If ^i, . . . , ^5 are any five tetrahedra such that dih(5'i) + • • • + 
dih(55) > 271, then 

r(5i) + --- + r(55) <4.52pt 

This is a consequence of two other calculations. 

Calculation 9.6.1. F < -0.207045 dih +0.31023815, provided dih e [do, 27r-4do], 
where do = dih(^(2, 2, 2, 2, 2, 2)). 

Calculation 9.6.2. F < 0.028792018, if dih > 27r - 4do. 

Calculations 9.1, 9.6.1, and 9.6.2 give a piecewise linear bound £ on F as a function 
of dihedral angle. See Diagram 9.6.3. (Again, there are minute discontinuities 
that may be eliminated in the same manner as before.) We claim that A.b2pt > 

£{ti)-\ 1-^(^5) whenever ti-\ hts > 27r. As in Lemma 9.5, we may assume that 

ti > do. Since £ is decreasing on [cZq, dihmax]) we may assume that + - • • + t5 = 27r. 
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Thus, only the interval [do, 2tt — Ado] is relevant for the optimization. On this 
interval, the bound is linear, so i{ti) H h £(^5) < 5£{27t/5) < 4.52 pt. 




Calculation 9.7. T < 0.389195 dih -0.435643, if yi > 2.05. 



The next inequalities relate the solid angles to the compression. 
Calculation 9.8. T < -0.37642101 sol +0.287389. 
Calculation 9.9. T < 0.446634 sol -0.190249. 

Calculation 9.10. T < -0.419351 sol +0.2856354 + 0.001, if e [2,2.1], for 
i = 1,... ,6. 

The following calculation involves the circumradius. We leave it to the reader to 
check that the dimension-reduction techniques of Lemma 8.7.1 may still be applied. 

Calculation 9.11. F < -0.419351 sol +0.2856354, provided that yA,y5,yQ > 2.1 
and that the circumradius of S is at most 1.41. 

Calculation 9.12. T < -0.419351 sol +0.2856354 - 5(0.001), if > 2.1 for some 
z, and y4 e [2, 2.1]. 

Calculation 9.13. T < -0.65557 sol +0.418, if e [2, 2.1], for z = 1, . . . , 6. 
Calculation 9.14. sol(S') > 0.21. 



Calculation 9.15. T + {K - sol)/3 < 0.564978 dih -0.614725, where K = {An - 



c /lo^/10 ;J — ] „. ^ To oncl ] „. 



To o ol 
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Calculation 9.16. dih > 0.98 and sol > 0.45, provided yi G [2, 2.05], and 2/2, Us e 
[2,2.2]. 

Let r(S') be replaced by vor(5') in each of the Calculations 9.1 - 9.16 to obtain a new 
list of inequalities 9.1'- 9.16'. (In 9.11', we drop the constraint on the circumradius 
of S.) We claim that all of the inequalities 9.*' hold whenever 5" is a quasi-regular 
tetrahedron of circumradius at least 1.41. In fact. Inequalities 9.1', 9.2', 9.4', 9.5.1', 
9.5.2', 9.5.3', 9.6.1', 9.6.2', and 9.7' follow directly from 9.17 and the inequalities 
dihmin < dih < dihmax- Calculations 9.3, 9.14, and 9.16 are independent of F, and 
so do not require modification. Inequalities 9.8', 9.11', and 9.12' also rely on 9.17, 
9.18 and 9.19, inequality 9.9' on 9.14, and inequality 9.15' on 9.16. Inequalities 
9.10' and 9.13' are vacuous by the comments of 8.2.4. 

Write S = S {1/1,1/2, 1/3, 1/4,, 1/5, ye), vor = vor(S'), and let rad = rad(S') be the cir- 
cumradius of S. 

Lemma 9.17. If the circumradius is at least 1.41, then vor < — 1.8 pt. 

Proof: If sol > 0.91882, then the lemma is a consequence of Lemma 9.18. (The 
proof of Lemma 9.18, under the restriction sol > 0.91882, is independent of the 
proof of this lemma.) Assume that sol < 0.91882. 

If S' and S are Delaunay stars, related as in Section 8.7, then vor(S") > vor(S') 
because Sq is obtained by slicing a slab from So- This means that we may apply 
the dimension-reduction techniques described at the beginning of this section, un- 
less the deformation decreases the circumradius to 1.41. An interval calculation 
establishes the result when + 2/2 + 2/3 > 6.3 and the circumradius constraint 
is met (Calculation 9.17.1). When the dimension reduction techniques apply, the 
dimension of the search space may be reduced to four, and an interval calculation 
similar to the others in this section gives the result (Calculation 9.17.2). 

We lacked the computer resources to perform the interval analysis directly when the 
dimension-reduction techniques fail and found it necessary to break the problem 
up into smaller pieces when |/i + 1/2 + 2/3 < 6.3. We will make use of the following 
calculations. 

Calculation 9.17.1. vor(S') < -1.8 pi provided 2/1 + 2/2 + 2/3 > 6.3, rad(5') = 1.41, 
and so\{S) < 0.91882. 

Calculation 9.17.2. vor(S') < —1.8 pt provided that rad(S') > 1.41, 2/1+2/2 + 2/3 > 
6.3, so^S") < 0.912882, and S lies in one of the seven subspaces of smaller dimension 
associated with /(sol, F) at the beginning of the section. 
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Calculation 9.17.3.1. Assume that yi +2/2 + I/3 < 6.3 and that rad > 1.41. Then 
sol > 0.767. 

Calculation 9.17.3.2. Assume that j/i + j/2 + J/3 < 6.192 and that rad > 1.41. 
Then sol > 0.83. 

Calculation 9.17.3.3. Assume that 2/1 + 2/2+2/3 < 6.106 and that rad > 1.41. 
Then sol > 0.87. 

Calculation 9.17.3.4. Assume that 2/1 + 2/2 + 2/3 < 6.064 and that rad > 1.41. 
Then sol > 0.9. 

Calculation 9.17.3.5. Assume that 2/1 + 2/2+2/3 ^ 6.032 and that rad > 1.41. 
Then sol > 0.91882. 

In the interval arithmetic verification of Calculations 9.17.3, we may assume that 
rad = 1.41 and that 2/1 + 2/2 + 2/3 is equal to the given upper bound 6.3, 6.192, etc. 
To see this, we note that the circumradius constraint is preserved by a deformation 
of S that increases 2/1, 2/2, or 2/3 while keeping fixed the spherical triangle on the 
unit sphere at the origin cut out by S. We increase 2/1, 2/2, and 2/3 in this way until 
the sum equals the given upper bound. Then fixing 2/1, 2/2, 2/3; and one of 2/4, 2/5; 
and 2/6, we decrease the other two edges in such a way as to decrease the solid angle 
and circumradius until rad = 1.41. 

This deformation argument would break down if we encountered a configuration in 
which two of 2/4, 2/5; and 2/6 equal 2, but this cannot happen when rad(5') > 1.41 
because this constraint on the edges would lead to the contradiction 

1.41 < rad{S) < rad(^(2.3, 2.3, 2.3, 2, 2, 2.51)) < 1.39. 



It is possible to reduce Calculations 9.17.3 further to the four-dimensional situation 
where two of 2/4, 2/5 and 2/6 are either 2 or 2.51. Consider a simplex S with vertices 
0, vi, V2, and vs. Let pi be the corresponding vertices of the spherical triangle cut 
out by S on the unit sphere at the origin. Fix the origin, V2, and V3, and vary 
the vertex v^. The locus on the unit sphere described as pi traces out spherical 
triangles of fixed area is an arc of a Lexell circle C. Define the "interior" of C to be 
the points on the side of C corresponding to spherical triangles of smaller area. The 
locus traced by vi on the circumsphere (of S) with \vi \ constant is a circle. Let C, 
also a circle, be the radial projection of this locus to the unit sphere. Define the 
"interior" of C be the points coming from larger \vi\. The two circles C and C meet 
at pi , either tangentially or transversely. The interior of C cannot be contained in 
the interior of C because the Lexell arc contains P2, the point on the unit sphere 
antipodal to p2 [FT, p. 23] , but the interior of of C does not. Furthermore, if |f 1 1 =2 
and 1^2! + |'i'3| > 4, then V2 or lies in the interior of C and C, so that the circles 
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have interior points in common. This means that one can always move vi in such a 
way that the sohd angle is decreasing, the circumradius is constant, and the length 
1^1 1 is decreasing (or constant if |vi| =2). If any two of j/4, j/5, and i/q are not at an 
extreme point, this argument can be applied to vi, V2, or V3 to decrease the solid 
angle. This proves the reduction. 

We are now in a position to prove the lemma for simplices satisfying yi + y2+y3 < 
6.3. Calculation 9.17.3.1 allows us to assume sol > 0.767. As in Section 8.6, let 
Sy — S{2, 2, 2, y, y, y). We will rely on the fact that YOT^Sy) is decreasing in 2/, for 
y e [2.26, 2.41]. In fact, the results of Section 8.6 specialize to the formula 

(9.17.4) vor 6„ = ^ttttttt ^rr H — arctan , 

^ ' ^ ^' (12-y2)i/2(i6_y2) 3 V 64-6y2 J' 

and the sign of its derivative is determined by a routine Mathematica calculation. 
It is clear that sol{Sy) is continuous and increasing in y. Since sol(S'2.26) < 0.767 
and sol(5'2.4i) > 0.91882, our conditions imply that sol(S') = sol{Sy) for some 
y e [2.26,2.41]. Let r(a) and vol(i?(a, b, c)) be the functions introduced in Section 
8.6.4. 

This suggests the following procedure. Pick y so that sol(S') > so^^Sj^). Calculate 
the smallest (or at least a reasonably small) a for which 

e(z/, a) := vov{Sy) - A5oct{l - l/a^)2r{a) 

is less than — 1.8 pt. Monotonicity (Calculation 9.20.1) and Lemma 8.6.5 imply 
that vor(5') < — 1.8 pt, if 2/1+2/2 + 2/3 > 2(2 + a). To treat the case that remains 
(2/1 + 2/2 + 2/3 < 2(2 + a)), use Calculations 9.17.3 to obtain a new lower bound on 
so^/S), and hence a new value for y. The procedure is repeated until a = 0.016. 
Calculation 9.17.3.5 completes the argument by covering the case 2/1 + 2/2+2/3 < 
6.032. We leave it to the reader to check that 

^(2.2626,1.096), ^(2.326,1.053), ^(2.364,1.032), ^(2.391,1.016) 

are less than — 1.8 pt and that 

SOl(>52.2626) < 0.767, SOl(<S2.326) < 0.83, SOl(52.364) < 0.87, SOl(<S2.39l) < 0.9. 

This completes the proof of Lemma 9.17. I 



Lemma 9.18. If rad(,S) > 1.41, then vor < -0.419351 sol +0.2856354. 

Proof: We adopt the notation and techniques of Lemma 9.17. If sol < 0.918819, 
then the result follows from Lemma 9.17. (The proof of 9.17 is independent of 

— i- c^^^ ] — — — — i;j ,1 — ^ t „j- i"/ c»^ 
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-0.419351 sol(^) + 0.2856354 - vor(^). We show that f{S) is positive. We will 
use the inequality f{S) > /(<S'tan) + 45octVol(S'\S'tan) of Lemma 8.6.4. A routine 
calculation based on the formula 9.17.4 shows that f{Sy) is increasing, for y e 
[2.4085,2.51]. 

If sol > 0.951385, then we appeal to Fejes Toth's convexity argument described 
in 8.6. To justify the use of his argument, we must verify that the cone over S'tan 
contains the circumcenter of S. The first, second, and third edges have length 2, 
and the fourth, fifth, and sixth edges are between 2.21 and 2.51 (Calculation 9.18.2). 
These are stronger restrictions on the edges than in the proof of Lemma 8.6.5, so 
the justification there applies here as well. We observe that sol(S'2.4366) < 0.951385 
so that 

f{S) > /(-Stan) > f{Sy) > /(-S2.4366) ~ 0.000024 > 0, 

where y satisfies sol{Sy) = sol(5'tan)- We may now assume that 0.918819 < sol < 
0.951385. 



Lemma 9.18.1. The combined volume of the two Rogers simplices along a common 
edge of a quasi-regular tetrahedron S is at least 0.01. 

Proof: The combined volume is at least that of the right-circular cone of height 
a and base a wedge of radius y/b'^ — a? and dihedral angle dihmin) where a is the 
half-length of an edge and 6 is a lower bound on the circumradius of a face with an 
edge 2a. This gives the lower bound of 

|(fe2 _ a2)a^^ > 0.1439(62 - a"). 

We minimize h by setting 6^ = ?7(2,2,2a)2 = 4/(4 — o?). Then 6^ — = (a^ — 
2)2/ (4 — a^), which is decreasing in a e [1, 2.51/2]. Thus, we obtain a lower bound 
on iP' — by setting a = 2.51/2, and this gives the estimate of the lemma. I 



We remark that sol (5*2. 4085) < 0.918819. If 1.15 < (2/1 2/2 +2/3 -4) /2, then Lemma 
9.18.1 and Section 8.6 give 

fiS) > /(-52.4085) + 45oct(l - 1/1.15^)0.01 > 0. 

(Analytic continuation is not required here, because of the constraints on the edges 
in Calculation 9.18.2.) We may now assume that 2/1+2/2+2/3 < 6.3. To continue, 
we need a few more calculations. 

Calculation 9.18.2. If sol > 0.918, then ^4,^5, 2/6 > 2.21. 

Calculation 9.18.3.1. If yi + 2/2 + ys < 6.02 and rad > 1.41, then sol > 0.928. 
Calculation 9.18.3.2. If ui + + < 6.0084 and rad > 1.41. then sol > 0.933. 
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Calculation 9.18.3.3. If yi+y2 + y^< 6.00644 and rad > 1.41, then sol > 0.942. 

In the verification of these calculations, we make the same reductions as in the 
Calculations 9.17.3. 

Adapting the procedure of Lemma 9.17 and Lemma 8.6.5, we find that a lower 
bound on the solid angle leads to an estimate of a constant a with the property 
that f{S) > whenever yi+y2 + y3 > 4 + 2a. That is, we pick a so that a) := 
f{Sy) + 85oct(l — l/a^)r{a) is positive, where y is chosen so that sol{Sy) is a lower 
bound on the solid angle. The values 

^1(2.4085,1.01), ^1(2.4165,1.0042), ^i (2.42086, 1.00322), ^i (2.4286, 1.0017) 

are all positive. This yields the bound 2/i + 2/2 + 2/3 < 6.0034. 

Assume that S satisfies 2/1+2/2 + 2/3 < 6.0034 and rad(5') > 1.41. Then by Cal- 
culation 9.18.3.3, sol(S) = so^^tan) > 0.942. Also rad(^tan) > rad(;5)/1.0017 > 
1.41/1.0017 > 1.4076, because rescaling S'tan by a factor of 1.0017 gives a sim- 
plex containing S. This means that S'tan satisfies the hypotheses of the following 
calculation. Calculation 9.18.4 completes the proof of Lemma 9.18. I 

Calculation 9.18.4. If rad(5tan) > 1-4076 and 0.933 < sol(5tan) < 0.951385, then 

/(^tan) > 0. 

In this verification, we may assume that the fourth fifth, and sixth edges of S'tan 
are at least 2.27. for otherwise the circumradius is at most 

rad(5(2, 2, 2, 2.27, 2.51, 2.51)) < 1.4076. 

In the verification of 9.18.4, we also rely on the fact that 

fi{x4,xs„X6) := /(S'(2,2,2, V^, V^)) 

is increasing in (0:4, 0:5, xg) G [2.27^,2.51^]. Here is a sketch justifying this fact. 
The details were carried out in Mathematica with high-precision arithmetic. The 
explicit formulas of Section 8.6 lead to an expression for df 1/8x4 as 

W{X4, X5, Xq){-X4 + X5+ Xq) 

(-16 + a;4)2A3/2 

where is a polynomial in X4, X5, and xq (with 13 terms). To show that W is 
positive, expand it in a Taylor polynomial about (2:4, x^, xq) = (2.27^, 2.27^, 2.27^) 
and check that the lower bound of Inequality 7.1 is positive. 



r^„\ n -I n jc ^ 
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For a < 6 < c, we have vol(i?(a, b, c)) = o((62 - a'^){c'^ - b'^))^/^/6. The final four 
calculations are particularly simple (to the extent that any of these calculations are 
simple), since they involve a single variable. They were verified in Mathematica 
with rational arithmetic. 

Calculation 9.20.1. The function (l-l/a^)vol(i?(a, r]{2, 2, 2a), 1.41)) is increasing 
on [1,1.15]. 

Calculation 9.20.2. The function vol(i?(a, 7/(2, 2, 2a), 1.41)) is decreasing for a e 
[1,1.15]. 

Calculation 9.20.3. For a e [1, 1.15] we have 

vol{R{a, r]{2, 2, 2a), 1.41)) < vol{R{a, r]{2, 2.51, 2a), 1.41)). 



Calculation 9.21. Let pi — pi{yi, . . . ,2/6) be the point in Euclidean space intro- 
duced in Section 8.6.7. For y e [2.3, 2.51], 

^(77(2.51, 2.51, y) - y'j^fl'' > -0.75 > ^{\pi{y, 2.51, 2, 2, 2, 2.51)1^ - y^/A^/^ 
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Appendix. Proof of Theorem 6.1 

Douglas J. Muder 

Notation and observations. Let P be a point of degree d. If we consider P to 
be the center of the configuration, then the first rim of points wiU be the set of d 
points adjacent to P, and the second rim will be those at distance 2 from P. Let 
5{P) be the sum of the degrees of the first rim points. If d = 6, it is easy to see 
that the number of points on the second rim is s = 5{P) — 24, and the total number 
of points of distance at most 2 from P is 5{P) — 17. The second rim is thus an 
s-gon. This s-gon and the triangulation of the P-side of the s-gon will be referred 
to as the inner graph. The inner-graph degree of a second rim point will be called 
its inner degree. The number of second-rim points with inner degree 5 is equal to 
the number of degree 4 points in the first rim, and the number of second-rim points 
with inner degree 3 is equal to the number of degree 6 points on the first rim. All 
other second-rim points have inner degree 4. Points beyond the second rim will be 
called extra points. If there are no extra points, then there are at least 2 second-rim 
points whose degrees are equal to their inner degrees. Let N^{P) be the number 
of points of degree A in rim i around P. Notice that if d = 6, then 

6{P) = 30 -Nl{P) + N^{P). 

Further, Eulcr's formula gives N — 12 = Nq — N4. Putting this together with 5.1.1 
and 5.1.5, we see 

13 < 12 + A^6 - A^4 < 15 
1 < 1 + iV4 < iVe < 3 + iV4 < 5. 

Lemma 1. There is no triangle whose vertices all have degree six. 

Lemma 2. Suppose Nq > 3. Let Gq be the subgraph of points of degree six. Either 
Ge is three points in two components, or it is a single (open or closed) path, with 
no other edges. 

With these lemmas and this notation, we consider the difi^erent possible values of 
iVe- 

iVe = 1. This forces N4 = 0. Let P be the degree 6 point. Then 6{P) = 30, aU 
13 points are in the inner graph, and all 6 second-rim points have inner degree 
4. But at least two of these second-rim points have inner degree equal to degree. 
Contradiction. 
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A^6 = 2. If the two degree 6 points are adjacent, let P be either one. Now 5(P) = 
31 — Nl, and we see that there are no extra points, and all second rim points have 
degree 5. All but A^^ second rim points have inner degree less than 5, and at least 
two of them will have inner degree equal to degree. So A4 > Nl = 2, and N < 12. 
Contradiction. 

If the degree six points are not adjacent, then the (6, 6) forces N > 14, so N4 — 0. 
If 5{P) is either of the two degree six points, 5{P) — 30, so the second rim has 6 
points, all of inner degree 4, and there is one extra point. If the extra point has 
degree 5, it is adjacent to all but one of the second rim points. This unique second- 
rim point is then part of a quadrilateral, which can only be triangulated by the 
diagonal edge that does not include the extra point. This creates two second- rim 
points of degree 6, in violation of our assumptions. Therefore the extra point has 
degree 6, and we have 6.2. 

Nq = 3. Either Gq is an open path or it has two components. In either case there 
is a (6, 6), so iV > 14 and iV4 = or 1. 

In the first case let P be the center of the path. Now 5{P) = 32 — Nl. We see that 
there are no extra points and any degree 4 point is in the first rim. The second 
rim has at most one point of inner degree greater than 4, and two points whose 
inner degrees are equal to their degrees. So there must be a point of degree 4 on 
the second rim. Contradiction. 

In the second case 5.1.8 forces N — 15, hence N4 — 0. If P is one of the points in 
the two-point component of Ge, then d{P) — 31, so there is one extra point and 7 
second-rim points. Either one or two of the second rim points is not connected to 
the extra point. In either case at least one of these two has degree = inner degree, 
which is at most 4. Contradiction. 

Nq = 4. Now Gq is either an open or closed path of length 4. If open, then by 
5.1.8, N = 15 and A4 = 1. If P is either of the interior degree 6 points, then 
d{P) = 32 — Nl, the second rim is an (8 — Nl)-gon with Nl points of inner degree 
at least 5 and Nl extra points. We cannot have Nl = 0, or else there would be at 
least two second rim points of degree 4, so Nl = 1. This holds no matter which 
interior degree 6 point we started with, so the degree 4 point must be adjacent to 
both. It also must be adjacent to one of the other two, or else there is a (6, 6,4). 
Therefore there is a degree 6 point whose first rim degree sequence is 6/4/6/5/5/5, 
producing a second rim inner degree sequence of 4/3/5/3/4/4/4. The point of inner 
degree 5 is the fourth neighbor of the degree four point, and cannot be degree 6 
without making Gq closed. Therefore there is an edge connecting the two points 
of inner degree 3. Adding this edge to the inner graph creates a hexagon with all 
points of "inner" degree 4, and the extra point in its interior. The extra point is 
connected to 5 of these Doints. and the sixth has deeree 4. Contradiction. 
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If Gq is a closed path, then a degree 4 point Q must triangulate this quadrilateral. 
Using Q as the center, there are 4 first rim points and 8 second rim points. The 
second rim can be drawn in a square with the 4 points of inner degree 3 on the 
corners and the 4 points of inner degree 4 at the midpoints. There is at most one 
additional degree 4 point (other than Q). If it is anywhere but at one of the second- 
rim midpoints, a (6, 6, 4) exists. The four midpoints cannot be adjacent to any 
additional second-rim points without creating an illegal triangle or quadrilateral. 
Since there can be at most one more degree 4 point, at least 3 of the 4 second- rim 
midpoints must be joined to the extra point(s). Two midpoints of consecutive sides 
can't be joined to the same extra point without forcing the corresponding corner to 
be degree 4, which is impossible. So there must be two extra points, and thus no 
degree four points other than Q. One extra point must be joined to the midpoints 
of each pair of parallel sides of the square. But all these edges cannot be drawn 
without intersecting. 

Nq = 5. Now Gq is a path of length 5, closed so that a (6, 6, 6) doesn't exist. Any 
point in this pentagon is nonadjacent to two others, so the situation of 5.1.8 applies 
in 5 different ways. Let P, Q, -R € Gq with P and Q adjacent, but neither adjacent 
to R. Then there are 8 points adjacent to either P or Q, and 6 points adjacent 
to R, so there are exactly two common points in these two sets. In our case these 
two points are precisely the other two points of degree 6. This means that no point 
inside or outside the pentagon can be connected to more than two points of the 
pentagon. Thus the only possible configuration for the 15 points is to form three 
concentric pentagons, with Gp, as the middle one. Prior to triangulating the inner 
and outer pentagons, all of their points have degree 4. Triangulating either creates 
more degree 6 points. Contradiction. 

Proof of Lemma 1. Let Pi, P2, and P3 be the vertices of such a triangle. If Pj 
has no neighbors of degree 4, then the only possibility is 5{Pj) — 32, which forces 
A^5 = 4, A/g = 2, and N = 15. All second-rim points have inner degree at most 
4, and there are no extra points. Therefore two second-rim points have degree 
and inner degree 4. Call them Qi and Q2- Therefore, Ar4 > 2, Nq > 5, and 
Nq > 2. Let Ri and R2 be second-rim points of degree 6. Each Qk and Rc must be 
adjacent, or else {Pj, R^, Qk) is a (6,6,4). But all edges from Qj are inner-graph 
edges. So Q1R1Q2R2 must be a second-rim quadrilateral, which is impossible. 
Thus Nl{Pj) > 1 for each Pj. 

This forces N4 = 2. At least one of the degree 4 points is adjacent to two of the 
Pj. One of these two, say Pi, is not adjacent to the other degree 4 point. Since 
5{P) < 32, we are left with four possible first-rim degree sequences for Pi: (1) 
6/6/4/5/5/5; (2) 6/6/4/6/5/5; (3) 6/6/4/5/6/5; or (4) 6/6/4/5/5/6. 

The last three possibilities are easily dealt with. Sequence (3) contains a (6,6,4). 
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vertices of the edge, and P3 and P4 be the other vertices of the triangles. The 
conditions Nl{Pj) > 1 and N4 = 2 force P3 and P4 to have sequence (1). In (2) 
we can invoke 5.1.8 with the 3 degree 6 points in the first rim. Call them R2 
and T, with Ri adjacent to R2- In order for the numbers to work out, there can 
be precisely two points which are adjacent to T as well as one of the Rj. But in 
(2) the center, the first rim degree 4 point, and the second rim neighbor of the first 
rim degree 4 point must all fit this description. 

Sequence (1) is more difficult to eliminate. Let Ri, . . .Re be the first rim points 
listed in the order of (1). Let Si, . . .,87 be the second rim points, with being 
adjacent to both Ri and R2, and having only Ri as a first rim neighbor. There 
is at least one degree 6 point in the second rim — if it is at ^i, then 5{Ri) = 32; if 
it is anywhere else we can invoke 5.1.8. In either case, = 15 and there are two 
degree 6 points outside the first rim. These must lie in { 5*6, 5*7, 5*1, 5*3} to avoid 
making a (6, 6, 4) with Ri and -R3. At most one of these points can be adjacent to 
i?i (since 5{Ri) < 33), so S3 must have degree 6. Now ^3 must connect to both the 
degree 6 and degree 4 points in {Sq, Sr, Si} to avoid making a (6, 6, 6) or (6, 6, 4) 
with the center point. But S3 has inner degree 5 and so can connect to at most 
one of those points. Contradiction. I 

Proof of Lemma 2. If Gq has at least three points and no triangles, then there 
exists a (6, 6) forcing N > 14. Any point in Gq can be adjacent to at most two 
others, or else either a triangle or a (6, 6, 6) is created. Therefore each component 
is a path. There cannot be three or more components without producing a (6, 6, 6). 
If Ge has two components, each of them must be a complete graph for the same 
reason. Since there are no triangles, both of the components must have at most 
two points. If there are two components with two points each, then 5.1.8 forces 
= 15. Each component is adjacent to 8 non-Ge points of the original graph, 
and not adjacent to 3. Now 5.1.8 forces each point in the other component to be 
adjacent to these 3. But two points can have only two common neighbors. So the 
only two-component Gg is the one described in the lemma. I 



